The Pan-STARRS Moving Object Processing System 

Larry Denneau^ (denneau@ifa.hawaii.edu), Robert Jedicke^, Tommy Grav^, Mikael 
Granvik^, Jeremy Kubica*^, Andrea Milani^, Peter Veres^, Richard Wainscoat^, Daniel 
Chang^, Francesco Pierfederici^, N. Kaiser^, K. C. Chambers^, J. N. Heasley^, Eugene. A. 
Magnier^, P. A. Price^, Jonathan Myers^, Jan Kleyna^, Henry Hsieh^, Davide 
Farnocchia^'^, Chris Waters^, W. H. Sweeney^, Denver Green^, Bryce Bohn^, W. S. 
Burgett^, J. S. Morgan^, John L. Tonry^, K. W. Hodapp^, Serge Chastel^, Steve Chesley^, 
Alan Fitzsimmons^^, Matthew Holman^^, Tim Spahr^^, David Tholen^, Gareth V. 
Williams^^, Shinsuke Abe^'^, J.D. Armstrong^, Terry H. Bressi^"^, Robert Holmes^^, Tim 
Lister Robert S. McMillan^^ Marco Micheli^ Eileen V. Ryan^^, WiUiam H. Ryan^^ 

James V. Scotti^^ 



-2 - 

Received ; accepted 



57 Pages, 26 Figures, 13 Tables 



^Institute for Astronomy, University of Hawai'i, 2680 Woodlawn Dr., Honolulu, HI 96822 
^ Johns Hopkins University, Baltimore, MD 

^Department of Physics, P.O. Box 64, 00014 University of Helsinki, Finland 
^Google, Inc. 

^University of Pisa, Pisa, Italy 

^Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 21218 
^Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 
^University of Arizona, Tucson, AZ 

^Jct Propulsion Laboratory, Pasadena, CA 

^^Astrophysics Research Centre, School of Mathematics and Physics, Queen's University 

Belfast, Belfast, BT7 INN, UK 

^^Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 

^^Smithsonian Astrophysical Observatory, Cambridge, MA 

^^Institute of Astronomy, National Central University, Taiwan 

^^Lunar & Planetary Laboratory, University of Arizona, Tucson, AZ 

^5 Astronomical Research Institute, 7644 NCR 1800E, Charleston, IL 61920 

■^^Las Cumbres Observatory Global Telescope Network, Inc., 6740 Cortona Dr. Suite 102, 

Santa Barbara, CA 93117 

^■^Magdalena Ridge Observatory, New Mexico Tech, 801 Leroy PL, Socorro, NM 87801 



-3- 



ABSTRACT 

We describe the Pan-STARRS Moving Object Processing System (MOPS), a 
modern software package that produces automatic asteroid discoveries and identi- 
fications from catalogs of transient detections from next-generation astronomical 
survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from 
a synthetic but realistic population of asteroids whose measurements were sim- 
ulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical 
grid population, we demonstrate that MOPS can detect populations of currently 
unknown objects such as interstellar asteroids. 

MOPS has been adapted successfully to the prototype Pan-STARRSl tele- 
scope despite differences in expected false detection rates, fill-factor loss and 
relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 
telescope and survey. MOPS remains highly efficient at detecting objects but 
drops to 80% efficiency at producing orbits. This loss is primarily due to con- 
figurable MOPS processing limits that are not yet tuned for the Pan-STARRSl 
mission. 

The core MOPS software package is the product of more than 15 person-years 
of software development and incorporates countless additional years of effort in 
third-party software to perform lower-level functions such as spatial searching or 
orbit determination. We describe the high-level design of MOPS and essential 
subcomponents, the suitabihty of MOPS for other survey programs, and suggest 
a road map for future MOPS development. 

Subject headings: Surveys:Pan-STARRS; Near-Earth Objects; Asteroids 
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Introduction 



As with most scientific endeavors, the history of asteroid and comet studies depicts an 
exponential increase in the rate of discovery since the identification of Ceres by Piazzi more 
than 200 years ago. This work describes the next step in the evolution of asteroid surveys 
— an integrated, end-to-end moving object processing system (MOPS) for the Panoram ic 



Survey Telescope and Ra pid Response System (Pan-STARRS 



2004 



Kaiser et al. 



2002 



Kaiser 



Hodapp et al. 



2004j ). The system's prototype telescope (PS-1) employs a 1.4 gigapixel 



camera able to detect asteroids and comets faster than ever before. 

The first asteroids were discovered serendipitously and laboriously by eye until 
the dedicated photo graphic surveys of the 1950s an d 1960s like the Yerkes- McDonald 



( jKuiper et al. 



19581 ) and Palomar Leiden Surveys (Ivan Houten et al.lll970l ). The 



photographic surveys required a major effort due to the need for human 'bli nking' of the 



Alvarez et al. 



i mage s to identify the moving objects. The realization in the early 1980s by 
( ll980l ) that the extinction of the dinosaurs ~65 million years ago was precipitated by the 
impact of a large asteroid or comet with the Earth stimulated enhanced funding for wide 
field asteroid surveys. While those surveys identified many asteroids and comets that may 



eventually strike t he Earth, only one ob 



2008 TC3 {e.g. 



Jenniskens et al. 



2009 



ect has been discover ed that actually hit the Earth 



Boattini et al. 



2009h . 



The photographic sear ches leveraged decades o: 



surveying but Spacewatch ( IGehrels 



1991 



McMillan 



expe rience in wide field astronomical 



20071 ) spearheaded the first use of 



CCDs in asteroid surveys. In the beginning, the small CCDs of the time limited the success 
of this search program but the asteroid discovery rate improved dramatically when they 
obtained a high quantum efficiency 2Kx2K CCD with an ~30 arcmin field of view and 
adopted a 'drift scanning ' survey technique that eliminated the need for long readout 



times. At the same time. 



iRabinowitz 



(119911) developed the first automated moving object 
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detection program to identify asteroids and comets in Spacewatch's drift-sca n images tha t 



launc hed the con temporary generat ion of wide field surveys such as NEAT (IHelin et al 



19971), LONEOS flBowell et al. 



19951 ) ■ LINEAR f lStokes et al 



2000) and CSS ((Larson 



2003). 



These modern wide field aste roid surveys were prompted and subsequently funded by 



20081 ) with the goal of identifying 90% of Near 



the NASA Spaceguard Program (iHarris 

Earth Objectj^ (NEOs) larger than 1 km diameter before the end of 2008. All the surveys 
broadly employ similar techniques for asteroid and comet identification. Their wide field 
cameras have large pixel scales (typically >1 arcsec) and image a field 3-5 times within 
about an hour. Their moving object detection software identifies 'sources' in each image 
and then spatially correlates the detections between the images to identify and remove 
stationary objects. All the 'transient detections' are then searched for consistency with 
a single object moving linearly across the sky at a constant rate of motion between the 
exposures. The constant motion requirement allows the software to achieve per-detection 
signal-to-noise {S/N) levels of ~1.5 to 3a. Sets of 3-5 linked detections (a 'tracklet') 
representing candidate real moving objects are reported to the Minor Planet Center (MPC), 
though some groups review the detections by eye before submission to reduce the false 



detection rate. Overall, thes e groups wer e 



NEOs before 10 June 2008 ([Harris 



20081). 



wildly success: 



Mainzer et al 



ul, ide ntifying ~79% of the >1 km 



( 1201 ll ) report that the Spaceguard 



goal of discovering 90% of the >1 km diameter asteroids was actually met some time later 
but before 2011. 

The CCD surveys described above were designed to discover large, extremely hazardous 



^Near Earth Objects are asteroids or comets with perihelion distances of < 1.3 AU. Most of 



Gladman et al. 



these objects are in unstable orbits with dynamical lifetimes of <10 Myrs {e.g. 
20001 ) and will quickly be ejected from the Solar System or impact the Sun or Jupiter but a 
small fraction will eventually strike Earth. 
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NEOs, not to characterize the size-frequency distribution (SFD) of various solar system 
populations. While the process of searching for the NEOs, the y discovered several objects 
interesting in their own right {e.g. Comet Shoemaker-Levy 9, 
fast moving objects, typically very small and nearby asteroids, 



Shoemaker ( 



1995 



Rabinowitz et al. 



Jenniskens et al. 



; very 



(119931 ): 



( 120091 )): and post facto 



the Centaur object (5445) Pholus; 2008 TC3, 
determinations of their surveying efficiency allowed new estimations of the SFD of many 
small body populatio ns {e.g. obiects with orbits entirely inte rior to Earth's orbit, l EOs 



Zavodnv et al. 



Centaurs, 



(|2008[): NEOs. 



Jedicke and Herron 



Distant objects. 



Larsen et al. 



Rabinowitzl (11993f ): Main Belt, 



Jedicke and Metcalfe 



19971) : Trans- Neptunian objects, TNOs, 



Larsen et al 



9981): 



mm-, 



(I2OO7I )). Despite the success of these programs there is room 



for improvement in surveying techniques, moving object detection software, and reporting 
and follow-up methods: 



Survey Technique and/or depth 

The surveys in operation in 2003 were not capable of meeting th e Spaceguard 



Goal due to their sky-plane coverage and/or limiting magnitude (jJedicke et al. 



20031 ) . These surveys concentrated on the region of sky near opposition because 



asteroids are generally brightest in that direction but objects on orbits that will 
impact the Earth are under-represented towards opposition. Their sky-plane density 
increases at relatiyely small ecliptic latitu des and solar elongations < 90deg {e.g. 



Chesley and Spahr 



2004 



Veres et al. 



2OO9I ). However, asteroids in this region of the 



sky are notoriously difficult to observe because of the limited time that they are 
above the horizon, the high air mass when they are visible, their large phase angle 
(reduced illumination of the visible surface) and faintness. Furthermore, the 3-5 
repeat visits/night to the same field is wasteful of survey time. More sky could be 
covered or the same sky could be imaged to greater depth if there were fewer visits to 
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the same field each night. 



Follow-up 

To maximize survey coverage and their discovery rate most of the surveys visit a 
field on only one night per lunation to obtain a set of 3-5 linked detections. The 
detections provide the sky plane location and velocity vector for the object rather 



motion {e.g. 



Rabinowitz 


1991; 


Jedicke 


1996) 



Os, the surve ys flag objects with unusual rates of 



19961 ) as NEO candidates and then reacquire 



the objects in special follow-up efforts that reduce the time available for discovery of 
new objects. Alternatively, they report the detections to the MFC who then perform 
a more sophisticated probabilistic analysis to determine the likelihood that the object 
may be a NEO. If the object meets their likelihood cutoff, it is posted on their 
NEO confirmation webpagcl and hopefully recovered by professional and amateur 
astronomers around the world. Given that telescope time is a valuable commodity 
and that the next generation of professional asteroid survey telescopes will discover 
more and even fainter objects, it is clear that the survey telescopes must provide their 
own target follow-up and it should be incorporated directly into their survey pattern. 

Moving Object Detection and Orbit Determination Software 

With limited follow-up resources and inefficiency in the ability to determine if a 
candidate asteroid is likely a NEO based on its magnitude and apparent rate of motion 
(its 'digest' score), it would be useful for a survey's moving object detection system 
to determine orbits from observations over multiple nights instead of merely reporting 
sets of detections (tracklets) to the MFC. This allows for improved discovery rates of 
NEOs that are 'hidden in plain sight', indistinguishable from main-belt asteroids by 



^http:/ /www. minorplanetcenter.net/iau/NEO/ToConfirm. html 
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their brightness and velocity vectors alone, and therefore scoring too low using the 
MFC's NEO digest scoring to warrant follow-up. 

Reporting 

The system of reporting sets of detections to the MFC has been very effective and 
the follow-up response of the international community of amateur and professional 
observers has been fan tastic as amply demonstrated by the case of 2008 TC3 



(jJenniskens et al. 



20091 ) ■ Still, the process of submission by email to the MFC, 
posting on the NEO confirmation page, downloading and acquisition by observers is 
unwieldy. Since there is little coordination between observers there is the potential 
for wasted time with multiple and unnecessary follow-up. A modernized, automated 
reporting system tuned specifically to the capabilities and interests of each follow-up 
observatory would be useful to coordinate the follow-up effort. 

» Measured System Efficiency 

Existing asteroid detection software packages do not monitor their own detection 
efficiency and accuracy. Instead, analysis of the survey data must be performed post 
facto. It would be convenient if the packages incorporated an intrinsic near real-time 
measure of both their efficiency and accuracy for the purpose of system monitoring 
and subsequent data reduction. 

Some of the problems itemized above have been addressed by other s urveys, in 



Gladman et al. 



( I2OO9I ) using software 



particular the well-charac t erized main belt program of 

developed by iFetit et al.l (]2004l ). They generate synthetic detections in the same images 



used to identify their TNO candidates to measure the detection efficiency and accuracy. 

The future of wide field moving object processing systems is embodied in software 
such as the Fanoramic Survey Telescope and Rapid Response System (Fan-STARRS; 
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Kaiser et al.ll2002l : iKaiserl |2004| ) Moving Object Processing System (MOPS), described in 
the remainder of this paper. MOPS incorporates state-of-the-art spatial searching, orbit 
computation and database management into a cohesive package. MOPS resolves many of 



the issues described above and the system is capable of very high efficiency and accuracy. 



Pan-STARRS and Pan-STARRSl 



The University of Hawai'i's Pan-STARRS project was formed in 2002 to design 
and build a next-g eneration distributed-aperture survey system called Pan-STARRS4 



( iKaiser et al. 



2OO2I ). The Pan-STARRS4 design incorporates four distinct telescopes on a 
common mount on the summit of Mauna Kea. Shortly after inception, the project began 
development and construction of a single- telescope prototype system called Pan-STARRSl 



( IHodapp et all 12004 : iMorgan et al. 



20061 ) on Haleakala, Maui to validate essential 



components such as optics, camera and software pipelines. As of 2013, the Pan-STARRS 
project is constructing a second telescope called Pan-STARRS2 essentially identical to 
Pan-STARRSl, at the same site on Haleakala. The combination of Pan-STARRSl and 
Pan-STARRS2 is called PSl-l-2, and it is expected that the two telescopes will be operated 
together as part of a single science mission. 

The Pan-STARRSl telescope on Haleakala, Maui, began surveying for asteroids in 
the spring of 2010. This 1.8 m diameter telescope has a ~7 deg^ field of view and a 



1997 



2004h with 



~1.4 gigapixel orthogonal transfer array (OTA) CCD camera (ITonry et al.l 
0.26" pixels. Pan-STARRSl's large aperture, field of view and <13 s between exposures 
allows imaging of the entire night sky visible from Hawai'i to r ~ 21.2 in about 5 or 6 
nights. In practice, 2 to 8 images of the same field are acquired each night and consequently 
it takes longer to cover the entire night sky. 
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The Pan-STARRSl telesco pe is operated by the Pan-STARRSl Science Consortium 
(PSISC; 



Chambers 



2006 



20071 ) with a nominal mission of 3.5 surveying years. The PSISC 
survey plan incorporates about a half dozen sub-surveys but only three are suitable for 
moving object discovery (all the surveys are suitable for moving object detection because 
they all obtain at least one pair of images of each field each night) — the Svr all-sky survey, 
the Medium Deep (MD) survey and the solar system survey. 

The Pan-STARRSl detector system employs six passbands: gpi, rpi, ipi, zpi, ypi and 
wpi. iTonry et al.l (l2012i ) provide a detailed description of the Pan-STARRSl photometric 
system — the first four passb ands were designed to have similar characteristics to the 



SDSS {e.g. 



Karaali et al. 



20051 ). ypi was designed to take advantage of the good sensitivity 
of the Pan-STARRSl camera near l/xm and the wpi has a wide bandpass with high 
and low wavelength cutoffs designed to optimize the S/N for S and C class asteroids. 
Transformations from each filter to the V band are provided in table [1] for an object 
with solar colors and also for an object with an average S-l-C class spectrum. The 
difference between the solar and mean asteroid class transformations are significant since 
the Pan-STARRSl p hotometric system currently provides better than 1% photometry 



( Schlafiy et al. 



2OI2I ) in at least the Sir survey's filters. 



The Sir survey images the entire sky north of —30° declination multiple times in three 
passbands {gpi, rpi, ipi) in a single year. About half the sky within ±30° (two hours) of 
opposition in R.A. is imaged in the different passbands each lunation. The solar system 
survey is performed in the wide filter (wpi) and most of that time is spent near the ecliptic 



and near opposition to maximize the NEO discovery rate. A small percentage of solar 



system survey time is used to image the 'sweet spots' for Potentially Hazardous Asteroidi 



idB 



^PHAs are NEOs with orbits that approach to within 0.05 AU of the Earth's orbit and 
have absolute magnitudes if < 22. 
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(PHAs) within about ±10° of the echptic and as close to the Sun as possible subject to 
altitude, sky brightness constraints, etc. The MD survey obtains 8-exposure sequences at 10 
different fixed footprints on the sky, with an integration time of 1920 seconds per sequence. 
MD sequences are observed using either 8 x 240-second exposures in a single filter {ipi, 
zpi, ypi), or back-to-back sequences of 8 x 120-second exposures in gpi and rpi. 

Visits to the same footprint within a night are usually separated by about 15 minutes, 
a 'transient tir ne interval' ("^ ^TI), suitable for asteroid detection. The Image Processing 
Pipeline (IPP; lMagnierll2006l ) produces a source list of transient detections, i.e. new sources 
at some position in the image and/or those that change brightness but are otherwise not 
identifiable as false transients like cosmic rays or other image artifacts. The IPP then 
publishes the catalogs of transient detections to MOPS which searches for moving objects. 



3. The Pan-STARRS Moving Object Processing System (MOPS) 



3.1. Overview 



MOPS software development began in 2003 under management by the Pan-STARRS 
Project, with a mission to discover hazardous NEOs and other solar system objects while 
providing real-time characterization of its performance. Initial funds for Pan-STARRS 
construction and engineering were obtained by the University of Hawai'i via a grant from 
the United States Air Force Research Laboratory (AFRL). MOPS was the first and only 
science client to be funded by the Pan-STARRS Project — systems engineering, testing, 
performance tuning and creation of the Synthetic Solar System Model (S3M, §3.101) 
were directly supported by project resources. Additional in-kind software contributio ns 



from external collaborato rs in the areas of o rbit determination (iGranvik et al. 



Milani and Gronchi 



2009 



20101), spatial searching (IKubica et al.l 120071 ) and algorithm design 



- 12 - 



( iTyson and Ange]||200lf ) were incorporated via memoranda of understanding (MOUs). 
In 2008, the NASA Near Earth Objects Observations (NEOO) Program Office began its 
support of continued MOPS development and Pan-STARRSl operations. 

MOPS is the ffist integrated detector system that processes data from per-exposure 
transient detection source hsts through orbit determination, precovery and attribution. 
We sought to create a system that could independently measure our system's end-to-end 
throughput and effi ciency for the purpose of correcting for observational selection effects 



( jjedicke et al. 



2OO2I ). Our philosophy was that it is more important to accurately know 
the system efficiency than to be place all our effort in optimization. Figure [T] provides 
a high-level view of data flow into, through and from MOPS. Each of the MOPS data 
processing sub-steps is described in the following subsections. 

In the Pan-STARRS system design there is a functional separation between the image 
processing (IPP) and movin g obje ct processing (MOPS). Unlike the object detection 



algorithm of e.g. 



Petit et al. 



(120041 ) ■ the MOPS design mandates no integration with the 
image data. This was an intentional design decision based on the scale and organization 
of the Pan-STARRS project. The IPP is responsible for monitoring and reporting its 
detection efficiency and accuracy as is MOPS. In this way each Pan-STARRS subsystem 
can be developed and characterized independently. 

MOPS data processing consists of assembling groups of transient detections into 
progressively larger constructions until there are enough detections to produce a high-quality 
orbit believed to represent a real asteroid — a 'derived object'. Intra-night (same night) 
detection groupings are called 'tracklets' and inter-night (multi-night) groupings are called 
'tracks'. All tracks created by MOPS are evaluated by an orbit computation module 
after which a track is either 1) deemed to represent a real asteroid, or 2) rejected and 
its constituent tracklets returned to the pool of unassociated tracklets. When a derived 
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object is created MOPS uses the object's computed (derived) orbit to search for additional 
detections of the asteroid in MOPS data to further refine the orbit. 

MOPS expects its incoming transient detection stream to be organized by exposure, 
or field in MOPS parlance. Each field within MOPS is defined by its metadata e.g. 
(right ascension, declination [RA, dec]) boresight coordinates, exposure date and duration, 
filter. A MOPS detection by definition occurs in exactly one field and is defined by 
its position e.g. (RA, dec) in the field, observed magnitude, signal-to-noise (S/N), and 
associated uncertainties. 

3.2. Pipeline Design 

The MOPS pipeline operates by using a linear nightly processing model where data are 
ingested and processed in the order they are observed. Nightly data are ingested from a live 
transient detection stream and discrete processing stages are executed until all processing 
is completed for the entire night. Many of the MOPS processing stages build upon data 
structures created from prior nights. 

Out-of-order processing can be handled in a limited number of modes of operation as 
necessitated by the current Pan-STARRSl system. Enhancements to the MOPS pipeline 
to perform full processing on out-of-order data while preserving essential MOPS efficiency 
computations is still under development. The fundamental difficulty in out-of-order 
processing lies with the amount of data that needs to be recomputed when new observations 
are inserted in the middle of the temporal data stream. The existing MOPS design prefers 
nights to be added incrementally; insertion into the middle of the existing dataset essentially 
forces all subsequent nights to be reprocessed. 

Pipeline resource scheduling and management are handled by the Condor high- 
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throughput-computing software (IThain et aLll2005l ). Condor provides effective, flexible 



management of hardware resources that need to run the MOPS production pipeline 
simultaneously with test simulations or experimental processing of MOPS data. 

The MOPS pipeline is designed to be reliable in the event of cluster resource limitations 
and hardware failure {e.g. power outages, node failures). Working in tandem with Condor's 
process management infrastructure, the MOPS pipeline can be restarted easily at the 
proper point in the pipeline with all data structures intact. As an example, in 2011 the 
MOPS production MySQL database suffered a complete failure yet the MOPS pipeline was 
back online within hours after restoring the database from a backup. 

The discussion of each element of the MOPS pipeline is deferred to ^ where we 
follow the processing of detections from beginning to end of the pipeline and quantify the 
performance of each step using Pan-STARRSl data. 



3.3. Hardware 

The Pan-STARRSl MOPS runs on a modest cluster of standard Linux rack- mounted 
computers. MOPS makes no special demands on hardware so long as the cluster can 
keep up with incoming data. During early stages of MOPS design we based hardware 
requirements on estimates for transient detections stored and orbits computed per night of 
observing for a Pan-STARRS4-like mission, and scaled down the storage for Pan-STARRSl 
volumes. The current Pan-STARRSl processing hardware is capable of keeping up with 
a Pan-STARRS4-like data stream except for the storage of detections in the database. 
Table [2] lists the major hardware components employed by the Pan-STARRSl MOPS 
production processing cluster, and Fig. |2] shows the functional relationships between the 
MOPS hardware components. 
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To ensure against most of the types of disk failures we are likely to encounter, we use a 
multiple-parity redundant array of independent drives (RA1D6) for our database storage, 
meaning multiple parity bits for every byte stored on disk. This allows up to two drives to 
fail concurrently and still keep the disk array operational. 

Condor's workflow management automatically detects host-level failures and 
redistributes jobs appropriately. An entire cluster computation node can fail during a 
processing step and the MOPS pipeline will adjust seamlessly. 

MOPS pipeline administrative functions are handled by a 'fat' workstation-class 
computer that runs the main MOPS pipeline. Condor manager and MOPS web interface 
(described later). We distribute these functions between two identically configured hosts so 
that one can fail completely, its responsibilities then assumed by the other node at slightly 
degraded performance. During Pan-STARRSl commissioning and operations we have 
experienced complete database host failures, accidental deletion of a production MOPS 
database by a MOPS administrator, and numerous compute node failures, without loss of 
data or significant interruption of MOPS processing. 

3.4. Database 

Table [3] shows the storage required for MOPS database components under 
Pan-STARRSl and estimates for Pan-STARRS4. The catalog of high-significance 
detections and derived data products are stored in an industry-standard relational database 
to maximize interoperability with external analysis software and provide data mining 
capability. Data are organized into multiple row-column tables in which rows of one table 
may be related to rows in another table (see fig. [3]). For example, the MOPS FIELDS table 
contains exposure metadata describing the telescope pointing, exposure time, filter, etc. 
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for a single exposure. Each FIELD record is additionally assigned a unique identifier called 
FIELD_ID. A second table called DETECTIONS contains all transient detections ingested 
by MOPS and, to relate detections to a particular exposure, all detections from the 
exposure are assigned a FIELD_ID upon ingest that matches their exposure's FIELD_ID . 
Subsequent derived relations {e.g. groupings of detections into an asteroid 'identification' 
and computation of its derived orbit) are similarly maintained and are described in detail 
later. 

A relational database allows MOPS data to be manipulated, analyzed and queried by 
any external software that supports structured query language (SQL). Requests for both 
a small amount of data (such as a detection) or a large amount of data {e.g. all derived 
objects with q < 1.3) are made through a SQL query and the results are returned in a 
tabular representation. 

The MySQL relational database management system (RDBMS) was selected for MOPS 
because of its combination of performance and cost but other database vendors such as 
PostgreSQL and Oracle were considered and evaluated early in MOPS development. Any 
modern database system would work for MOPS as they all scale to billions of rows and 
have additional features that promote data integrity and system reliability. 



3.5. Low Significance Database 

MOPS was designed to perform searches for recovery observations using a lower- 
significance data archive called the 'low-significance dataset' (LSD). In normal MOPS 
processing a nominal confidence level, typically 5a, is used as a baseline for all processing. 
The Pan-STARRS4 camera and image processing are expected to deliver ~ 1500 5a 
detections per exposure. Below this confidence level the density of false associations 
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becomes so great that it overwhelms the data processing. However, the search for recovery 
or precovery observations of an object with a known orbit may constrain the predicted 
position and velocity enough to make it possible to identify observations in the LSD. 

Due to the much larger number of lower-significance sources in an astronomical image 
(there are 1,000 x more 3a than 5a detections) LSD detections are not stored in the MOPS 
relational database. Instead, they are stored in flat files on the local or network file system. 
The LSD archive is designed to be compact and efficiently searchable by exposure epoch, 
right ascension and declination. During only the attribution and precovery phases of 
MOPS processing, candidate 3a recovery observations are extracted for analysis from this 
dataset based on predicted positions, and successful recoveries are than 'promoted' into the 
high- confidence database. 



3.6. Data Exchange Standard 



Significant work related to early MOPS development and simulations led to creation 
of the Data Exchange Standard (DES) for Solar System Object Detections and OrbitC- 
The DES describes a file format for dissemination of observations and orbits of solar 
system objects and is used within MOPS and betwee n MOPS and add-on MO PS software 



2008 



20121). It provides 



developed by A. Milani and the OrbFit Consortium (iMilani et al.l 
mechanisms for reproducing detection<^tracklet and tracklet<S=^derived object relations that 
exist in the MOPS database, and allows for specification of statistical uncertainties for 



observational measures. The DES includes provisions for propagating 



^see Pan-STARRS document PSDC-530-004 by Milani et ai, Data Exchange Standard 
(2.03) for solar system object detections and orbits: A tool for Input/Output definition and 
control. 
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• orbit covariance and normal matrices, 

• detailed per-detection residuals, 

• ephemerides for real and synthetic asteroids, 

• detailed tracklet metrics beyond the scope of MOPS, and 

• radar and spacecraft observations. 

3.7. Web Interface 

Early in MOPS development we realized the utility of visual, web-based interrogation 
of the MOPS database, leading to the creation of the MOPS web interface. The MySQL 
RDBMS provides software 'hooks' into its internal code that allow web-friendly scripting 
languages to perform queries so that web-based user interfaces and reporting tools can be 
developed rapidly. 

The MOPS database user has ready visual representation of the MOPS processing 
stream, from a sky map of the nightly survey pattern, detection maps of each observed field 
indicating both synthetic and real detections, to the orbit distributions of derived objects, 
to the efficiency and accuracy of each of the processing steps. Pages for a single MOPS 
derived object show the object's linking history and the evolution of its orbit at each step. 
Fig. m shows the nightly sky map of an active MOPS database. 

3.8. Data Export 

The primary consumer of MOPS exported data products is the lAU Minor Planet 
Center (MPC) that maintains the authoritative catalog of observations and orbits for 
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minor planets, comets and natural satellites. As of early 2013, the MPC prefers email data 
submission because the Pan-STARRSl survey is oriented toward reporting single-night 
tracklets, not multiple-night derived objects. 

Under the Pan-STARRS4 processing model where MOPS produces derived objects 
with secure orbits on a regular basis, MOPS is capable of exporting monthly DES catalogs 
of detections, tracklets, orbits and indentification records for known and unknown derived 
objects. Pilot submissions to the MPC showed that this export and publication model 
works as long as the false derived object rate is suitably low {e.g. < 0.1%; Spahr, personal 
communication) . 

Early Pan-STARRSl processing revealed a false derived object rate too high for 
automatic submission to the MPC and an unsatisfactory NEO discovery rate so, to 
improve the short-term science output and satisfy funding agency requirements, the PSISC 
abandoned the derived object processing model in favor of the traditional tracklet reporting 
method employed by other asteroid surveys. Processing of derived objects will be restored 
for the final post-survey MOPS processing but the derived object throughput will be small 
due to the Pan-STARRSl survey strategy being optimized for single-night NEO discovery. 

Detections are distributed to the MPC via email after a human 'czar' visually verifies 
that candidates are real using the MOPS web-based interface. (NEO candidates consume 
most of the czar's time because of the high rate of false tracklets at high rates of motion.) 
In this display, a page lists nightly probable real asteroid tracklets with highlighted 
NEO candidate tracklets (see Fig. E]). Tracklets that belong to known numbered and 
multi-opposition o 



(IMilani et al. 



ajects are identified by a MOPS add-on module called known_server 
2OI2I ) and automatically submitted to the MPC. A MOPS administrative 
email account receives MPC confirmation email to verify the MOPS submissions. MOPS 
maintains its own table of submitted tracklets to the MPC for easy reproduction of 
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submitted detection details and to prevent multiple submissions of the same tracklets. 



3.9. Alert System 

The MOPS Alert System is a software agent that runs independently of standard 
MOPS processing and searches for events that require immediate notification or follow-up. 
The specific type of event is user-configurable using a small piece of Python code written 
by a MOPS scientist or engineer to query for interesting events. A typical example would 
be a fast-moving near-earth object (NEO) candidate tracklet that is unlikely to be observed 
again by the survey telescope and/or needs immediate follow-up so that it is not lost. 
In this case, the alert system identifies the candidate tracklet by its sky-plane rate of 
motion and the fact that the tracklet 's positions cannot be attributed to a known asteroid. 
Additional examples of alerts are tracklets of unknown objects with comet-like motion or 
extended morphology, or known objects with magnitude anomalies that indicate activity. 

The alert system allows alerts to be triggered from either a MOPS derived object or 



an individual t racklet. The a 
VOEvent (e.g. ISeaman et al 



ert p ayload itself is an extensible markup language (XML) 



2011 



a standardized data structure that can be consumed 
by automatic data-processing agents. The alert deployment infrastructure can be changed 
easily; the current MOPS implementation uses both a restricted internet email distribution 
list and a private Twitter feed to PSl Science Consortium partners. 

Alerts are organized using a publish-subscribe paradigm in which multiple alert 
definitions can be published to various 'channels' that may be subscribed to by multiple 
users. For example, a 'NEO channel' might publish alerts for near-earth object alerts 
(g < 1.3 AU) and potentially hazardous objects (NEOs with H < 22.5 and minimum orbital 
intersection distance < 0.05 AU). Such a channel could be subscribed to by NEO follow-up 
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organizations. Alerts currently in use by the PSISC include unusual lightcurves within a 
single night and detections of objects from asteroid families of interest {e.g. Hildas). 



3.10. Synthetic Solar System Model 



O ne major MOPS innovation is the incorporation of a synthetic solar system model 



(S3M; 



Gray et al. 



20111 ) that allows us to monitor MOPS development and measure its 
performance in real-time operations as quantified by the metrics defined in §3.11 It is 
important to introduce realistic detections of solar system objects into the processing system 
to ensure that it can handle real objects and to me asure a reali s tic de tection efficiency. The 



Gravetal 



(I2OIII ) and we only provide 



solar system model has been described in detail by 
a brief overview here. 

The S3M is a comprehensive flux-limited model of the major small body populations 
in the solar system that consists of objects ranging from those that orbit the Sun entirely 
interior to the Earth's orbit to those in the Oort cloud. It even includes interstellar objects 
passing through the solar system on hyperbolic orbits. The S3M contains a total of over 13 
million synthetic objects from 11 distinct small body populations with the only requirement 
being that they reach V < 24.5 during the time period from roughly 2005-2015 (the 
anticipated operational lifetime of the Pan-STARRS survey). The time period requirement 
only affects a sub-set of the populations in the distant solar system. 

The S3M has proven invaluable for developing and testing MOPS but suffers from one 
major limitation — it only tests the software system for known types of objects. Since 
we cannot anticipate the properties of an unknown population we supplemented the S3M 
with an artificial 575,000-object 'grid' population. Grid objects have a random and flat 
distribution in eccentricity (0 < e < 1), sin (inclination) (0° < z < 180°) and in the angular 
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elements (0° < Q,u, M < 360°). The semi-major axes of the objects were generated using 
uniform random distributions over six different ranges as provided in table ID The ranges 
and number of objects were selected so that the sky-plane density of objects in each range 
is roughly equal and to provide denser coverage in the semi- major axis ranges that have 
large numbers of known objects. The absolute sky-plane density of the grid population is 
50 per field, a factor of 20-30 less than the realistic S3M on the ecliptic and much higher 
than the S3M off the ecliptic. 

The absolute magnitudes of the objects in the grid population were generated with 
uniform distributions in each semi-major axis range but with lower limits that decrease 
{i.e. the objects are made larger) according to the square of the semi-major axis so that the 
objects are visible from Earth. 

A major benefit of testing MOPS with the grid population is that they can appear at 
the poles where the sky-plane density of real or S3M objects is small. §4.31 discusses the 
results of simulations with the grid population. 

3.11. MOPS Efficiency Concepts 

During ~ 4 calendar years of MOPS development we continuously benchmarked its 
performance using 3 metrics for each MOPS subcomponent. The metrics count the correct 
and incorrect associations of synthetic sources available at each MOPS processing step. The 
metrics are: 

• Efficiency 

The fraction of available associations that were correctly identified. 

• Accuracy 

The fraction of associations that are correct, e.g. consisting of detections from the 
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same synthetic object definition. 

• Quality 

The fraction of correct associations that meet or exceed a pre-defined quahty metric. 

To provide the detailed accounting required to measure system efficiency and accuracy, 
data structures for each association created in a MOPS processing phase are decorated with 
a label describing their disposition after processing. These labels describe the detections 
incorporated in the resulting data structure: 

• AVAILABLE - the data structure should have been created in the processing phase and 
is available to be assigned. 

• CLEAN - the data structure contains only synthetic detections or tracklets that belong 
to the same synthetic object, i.e. is a 'correct' data structure. This label might 
describe a tracklet that is created from two detections of the same synthetic object, 
or a derived object consisting of three tracklets that are each CLEAN and belong to the 
same synthetic object. 

• MIXED - the data structure contains only synthetic objects but the detections or 
tracklets are from different synthetic objects. 

• BAD - the data structure contains both synthetic and nonsynthetic detections. For 
tracklets this is a common and normal occurrence since we expect that synthetic 
detections will occasionally fall near real detections. 

• UNFOUND - an expected data structure for a synthetic object was not created. UNFOUND 

data structures are 'dummy' structures that represent an operation that should have 
happened. For example, if two sequential fields at the same boresight contain two 
detections of the same asteroid but a tracklet is not created, a dummy tracklet 
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containing tracklet metadata is created so that the event can be counted and 
characterized. 

• NONSYNTHETIC - the data structure contains 'real' detections from actual telescope 
data. NONSYNTHETIC data consists of true detections of asteroids and false detections 
from image artifacts. 

The pipeline must 'peek' at the input and output data structures before and after a 
MOPS processing step to assign synthetic labels but this is the only time that the pipeline 
knows that a data structure is synthetic — the labels are never examined during any MOPS 
computation or algorithm. Data structures that are 'contaminated' — contain mixtures of 
real and synthetic detections or different synthetic objects — continue to be propagated 
through the pipeline. The only exception is for UNFOUND objects for which dummy structures 
are created and decorated with the UNFOUND label so that these occurrences can be counted 
by analysis tools. 

MOPS uses the labels to calculate the incremental efficiency after each processing 
stage. For a derived object, the state at each processing stage is preserved so that a 'paper 
trail' exists for each modification. This allows a complete stcp-by-stcp reconstruction of 
a derived object's history from its tracklets including the ability to indentify the exact 
processing step where an object was lost or incorrectly modified by MOPS. 

3.12. Simulations 

The MOPS pipehne is exercised through the creation and execution of a 'simulation'. 
Each simulation begins with a collection of telescope pointings called 'fields' and a 
population of synthetic solar system objects whose positions and magnitudes are computed 
for every field. Objects that 'appear' in fields are converted into detections and stored in 
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the MOPS database as though they were reported from an actual telescope. Detections 
can be 'fuzzed' astrometrically and photometrically using parameters that model the 
telescope performance; these include sky background, detector noise, and the point-spread 
function (PSF). The full fuzzing model incorporates a basehne astrometric uncertainty 
from the plate solutions added in quadrature with a flux-dependent positional uncertainty. 
Poisson-distributed (in S/N) false detections may be inserted in the fields to simulate noise 
from a real detector. After the injection of synthetic detections for a night the MOPS 
pipeline is invoked as though operating on real data. 

MOPS simulations are used both to verify correct operation of the pipeline and to 
interrogate performance of a hypothetical survey. A simple verification test consists of 
creating a small main belt object (MBO) simulation with several nights of synthetic fields 
spanning two lunations, executing the pipeline, and verifying that all expected MBOs were 
'discovered' through creation of a derived object and 'recovered' via attribution or recovery 
where possible. In this fashion all essential elements of the pipeline are exercised and 
verified after software modifications. 

For comprehensive testing, a larger, more realistic synthetic population of objects is 
inserted into the simulation, the MOPS pipeline is executed, and the output provides a 
quantitative assessment of the system's efficiency and accuracy. The pipeline is designed to 
be essentially 100% efficient at creating tracklets and derived objects for MBOs with the 
designed Pan-STARRS performance. 

When evaluating survey performance we configure MOPS to more precisely mimic 
the astrometric and photometric characteristics using filter- and field-dependent limiting 
magnitudes and FWHMs. Given the calculated V-band apparent magnitude for a synthetic 
asteroid, we calculate its expected signal-to- noise (S/N) using per-field detection efficiency 
parameters, and then generate a fuzzed magnitude and S/N. From the fuzzed S/N and the 
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field's FWHM we compute the asteroid's fuzzed riglit ascension and declination. If there is 
variation in sensitivity due to sky brightness or airmass, these effects can be modeled by 
supplying appropriate values for photometric zeropoint or sky noise in the field's detection 
efficiency parameters. 

During MOPS development we ran thousands of small- to medium-scale simulations for 
software unit testing and verification. More importantly, we ran several large simulations 
with synthetic solar system models containing millions of asteroids spanning several years 
of synthetic PSl observations. These simulations were run on our modest PSl MOPS 
production cluster and took many weeks or months to simulate the multi-year Pan-STARRS 
mission. For instance, our full_Slb simulation used the full S3M population in a 2 year 
simulation and required ~180 calendar days (see §4.ip while our four year NEO-only 
simulation (neo_4yr) required 88 calendar days (see §4.2p . At various times during the 
pipeline processing the cluster was offline for maintenance or power outages but MOPS is 
designed to be interrupted and restarted. 

The upstream image processing subsystem {e.g. the Pan-STARRS IPP) can improve 
MOPS fidelity by providing a Live Pixel Server (LPS) that reports whether a pixel at or 
near a specific (RA,dec) in an image falls on 'live' camera pixels. This service improves 
MOPS simulations in two ways: 1) by generating synthetic moving object data with 
behavior that simulates real observations {e.g. to account for the loss of field due to the 
camera fill factor) and 2) accounting for missing detections when attempting to recover 
known asteroids that are in the field of view. In a single Pan-STARRSl exposure there 
are many reasons an area of the detector can be inactive besides simple chip gaps, such as 
saturation from a bright object in the field of view, or a cell used for guide star tracking. 
Thus the LPS tells MOPS what areas of the detector are active for any given exposure. If a 
LPS is not available from the upstream image processing subsystem, MOPS has provisions 
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for allowing static specification of a detector mask that defines areas on a detector where 
detections cannot occur (see Fig. [H]) • Synthetic detections that land on a dead part of the 
detector (as determined by the LPS) are marked with a special 'unfound' decorator so that 
they can be accounted for but otherwise omitted from all subsequent processing. 

For increased simulation fidelity MOPS can accept a per-field detection efficiency for 
point sources and trailed detections as a function of S/N and rate of motion so that MOPS 
can generate more realistic synthetic photometry. 

4. MOPS Verification 

During MOPS development we undertook several large-scale simulations to probe the 
performance limits of the MOPS design and implementation and understand where further 
attention was needed in software development. These simulations consisted of a complete 
S3M or large S3M subpopulation, a multi-year set of telescope pointings, and detector 
performance that simulated a Pan-STARRS4 system. Our motivation was to investigate 
MOPS performance for the Pan-STARRS4 system as specified in its design requirements, 
not the reduced, changing performance of the single Pan-STARRSl prototype telescope 
undergoing development and commissioning. 

4.1. 2-year S3M simulation (full_Slb) 

The MOPS 2-year S3M simulation consists of the full S3M and 100,424 fields spanning 
245 distinct nights over nearly two years (29 Dec 2007 through 2009 Oct 22). The 
simulation took ~3 months to execute including various interruptions for cluster downtime 
and maintenance. It assumes a nominal positional measurement uncertainty of 0.01", a 
FWHM of 0.7", a fill factor of 100%, and a constant limiting magnitude oi R = 22.7, 



- 28 - 



representing a single idealized Pan-STARRSl telescope. Poisson-distributed false detections 
were added to each field at a density of ~200/deg^. The simulated field pointings were 
generated by the TAO survey schedulej^ configured to produce a Pan-STARRS-like survey 
covering ~3,600 deg^ near opposition and ~600 deg^ near each of the morning and evening 
sweetspots. Each area is visited three times per lunation, with some nights lost due to 
simulated poor weather. MOPS achieved an overall 99.99997% tracklet efficiency and 
99.26% derived object efficiency. Tables |5] and [6] detail the tracklet and derived object 
efficiencies per S3M sub-population. 

For classes other than NEOs, Jupiter-family comets(JFCs) and long-period comets 
(LPCs), MOPS has an efficiency better than 99% for converting three or more tracklets 
to a differentially corrected orbit when observed within a 14 day window. For NEOs 
our performance is not as good, but as stated elsewhere in this work, we have been 
more concerned with quantifying MOPS performance for various populations than 
in optimization. MOPS contains many runtime configuration parameters that allow 
performance of individual modules to be tuned, often at the cost of greater false results. 
Pan-STARRSl operations have led the project to a different mode of NEO discovery 
than originally envisioned, so we have not focused on improving NEO derived object 
performance. Similarly for LPCs, we can again improve efficiency by tuning MOPS orbit 
determination modules to reject fewer parabolic or hyperbolic orbits. 



Our differential correction performance ( §5.5p is probably overstated because the 
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Milani and Gronchi 



it initial orbit determination (lOD) software ( iMilani et al. 
20101) employed by MOPS can perform its own differential correction after lOD, and we 
elect to use it in this mode. Therefore orbits handed to the JPL differential correction 
module are almost always close to a minimum in the solution space where convergence 
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will occur. MOPS supports other packages besides OrbFit to perform lOD, and in early 
evaluation of these packages we saw outstanding performance from the JPL differential 
corrector. 

4.2. 4-year NEO simulation (neo_4yr) 

We ran the neo_4yr simulation to calculate upper limits on the NEO detection and 
discovery rates with a Pan-STARRSl-like system. This system differed from the full_Slb 
simulation by using a 0.1" baseline astrometric uncertainty and by containing almost 
four years of TAO simulated opposition and sweetspot observations over the expected 
Pan-STARRSl mission from March 2009 through January 2013. The simulation only used 
the 268,896 NEO subset of the S3M. We were not able to simulate the Pan-STARRSl 
focal plane fill factor, as the code had not been incorporated into MOPS yet, but we felt 
it would be reasonable to use a 100% fill factor and scale down our results based on true 
Pan-STARRSl detector performance. False detections were not added to this simulation to 
speed up execution time, but the full_Slb results showed that this would not impact our 
results. 

The neo_4yr simulation found over 9,405 derived objects of which nearly half had 
orbital arc lengths greater than 30 days. As with the 2-year full_Slb simulation the 
tracklet efficiency is essentially 100%, with 10 out of 105,439 tracklets lost. MOPS linking 
performance is similar to f ull_Slb (see table [7j) but we find that orbit determination 
performs less well, at 81.7%, purely due to greater RMS residuals across the orbit from the 
larger astrometric uncertainty. Again, this acceptance threshold can be relaxed to improve 
efficiency (fraction of correct linkages surviving orbit determination) at a cost of an increase 
in false linkages. 



- 30 - 



4.3. Grid Simulations 

We have run several grid simulations to ensure that we were not tailoring the efficiency 
to localized regions in (a, e, i) phase space. However, we do not inject the grid population 
into the detection pipeline during normal MOPS operations as this would needlessly 
increase the tracklet linking combinatorics in an artificial and unrepresentative manner. 

To illustrate MOPS operations over the entire phase space we ran the real 
Pan-STARRSl survey {i.e. actual pointings and cadence) for two consecutive lunations 
with the grid population and a realistic (0.1") astrometric error model. Figures [7] through 
[9] show that MOPS is ~ 100% efficient for semi-major axes in the range from 0.75 AU TO 
5,0000 AU, eccentricities in the range [0, 1] and all possible inclinations. Fifteen fast-moving 
objects were rejected (out of 190,973) by the rate limit of 5.0 deg/day in the current 
Pan-STARRSl production configuration to reduce false tracklets. 

Figures (TU] through [T2] show that the derived object efficiency is consistently high for 
all synthetic grid objects on orbits ranging from those with semi-major axes well within 
the Earth's orbit to beyond 100 AU, over all inclinations with no deterioration in the 
efficiency for retrograde orbits, and for nearly circular orbits to those with e < 1 . There is a 
degradation in the efficiency for e ~ 1 due to the orbit determination failing at the highest 
eccentricities. 

One important feature of the grid population is that it includes objects that can appear 
near the poles where the sky-plane density of real solar system objects is negligible, i.e. it 
allows us to check that MOPS works correctly even for objects at high declinations as shown 
in Figs. Uni and dH While tracklet creation efficiency is nearly 100% for all declinations 
< 88°, within 2° of the pole the grid simulation revealed that slightly conservative tracklet 
acceptance parameters rejected several fast-moving objects. The derived object efficiency 
is typically >95% for declinations ranging from —30° to the north celestial pole indicating 
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that the tracklet hnking algorithm works across the entire sky. The drops in efficiency to 
~90% are largely due to an inability to compute an initial orbit (see §5.6p which is in turn 
usually due to a sparse observation cadence for the track. 



5. MOPS Validation and Real Data 



Processing of data from the current generation of s urveys has been addressed by 



software pipelines already employed by Spacewatch {e.g. 



Rabinowitz 



1991 



Larsen et al 



20011 ). the Catalina Sky Survey, LINEAR, etc.. We did not want to reproduce their work. 
MOPS was designed for use with a high-quality transient detection stream delivered by a 
next-generation Pan-STARRS4 or LSST-class survey system with nearly 100% fill-factor 
and a relatively modest number of systematic false detections. In many respects, despite 
the much higher data volumes expected from these telescopes, moving object algorithm 
design is simpler due to the high quality astrometry and photometry of the input data. 
MOPS performs less effectively with real Pan-STARRSl data because it is the prototype for 
the next-generation Pan-STARRS4 survey — not a next-generation survey itself. 

That said, the Pan-STARRSl telescope represents a significant step forward in wide 
field survey astronomy, and has become a major contributor in NEO detection and discovery 
and the largest single-telescope source of asteroid observations while satisfying multiple, 
disparate survey programs. We have adapted MOPS to Pan-STARRSl to maximize asteroid 
detection and NEO discovery from the Pan-STARRSl data stream through the creation of 
tools to reject false detections and tracklets and allow for human review of data before they 
are submitted to the Minor Planet Center. Additionally, due to a realized survey cadence 
that is inefficient at producing enough repeat observations to compute reliable orbits, 
we have de-emphasized MOPS's derived object functionality in normal Pan-STARRSl 
processing. In practice, the MOPS pipeline is executed through the tracklet-creation stage. 
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with some "add-on" tools that are described below. MOPS is still used to its full capabihty 
as a simulation and research tool. 

In this section we describe in detail all the Pan-STARRSl MOPS processing steps and 
quantify its performance on the transients provided by the prototype telescope. 

5.1. Pan-STARRSl Surveying 

The Pan-STARRSl 'effective' camera fill factor (/) as measured by MOPS is ~75%. 
The loss of image plane coverage results from gaps between the OTAs, dead cells and gaps 
on the CCDs themselves, problematic pixels {e.g. anomalous dark current or non-linearity), 
and the allocation of sub-sections on the OTAs (known as cells) to high speed acquisition 
of bright guide stars. The original MOPS design assumed / ~ 100% (Pan-STARRS4) and 
three 2-detection tracklets/lunation for new object discovery so that with / < 100% the 
maximum achievable system efficiency (emax) for moving objects must be somewhere in 
the range of < tmax < f^, depending on the object's rate of motion, the distribution of 
'dead' image plane pixels/area and the telescope pointing repeatability. This implies that 
the maximum Pan-STARRSl MOPS end-to-end derived object efficiency for objects in the 
field of view must he in the range 26-75% even for objects that would be imaged at high 
S/N. 

In part to compensate for the fill-factor, in October 2010 the Pan-STARRSl survey 
was modified to obtain a four-exposure 'quad' in a manner similar to other NEO surveys 
instead of a single pair of exposures. The 4 images of a quad are typically acquired at the 
same boresight and rotator angle with a transient time interval (TTI) of about 15 minutes 
between sequential detections. Thus, a tracklet containing four detections would have a 
typical arc length of about 45 minutes. We then search for tracklets with three or four 
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detections in the quad. The false tracklet rate decreases as the number of detections in 
the tracklet increases but even with 4 detections per tracklet there are still false tracklets. 
To cope with the false tracklet rate we implemented a NEO 'czaring' procedure and an 
associated infrastructure for extracting 'postage stamp' 100 x 100 pixel images centered 
on every detection incorporated into a tracklet. A human observer, the 'czar', vets every 
tracklet before submission to the MPC. 

The tracklet processing for pair observations by other Pan-STARRSl surveys is 
specifically directed towards NEO detection and is divided into two rate of motion regimes 
— slow and fast, corresponding to 0.3 to 0.7 deg/day and 1.2 to 5 deg/day respectively. 
i.e. Tracklets containing just two detections are not detected if they have rates of motion 
between 0.7 and 1.2 deg/day. The reason for the two-regime processing is that the confusion 
limit from false detections becomes unmanageable around 0.7 deg/day, but beyond ~1.2 
deg/day we can use traihng information to reduce the false tracklet count. The slow 
regime's lower limit corresponds roughly to a rate of motion that easily distinguishes 
between NEOs and main belt asteroid motions at opposition. The upper limit in the slow 
regime is set empirically at the rate at which the confusion limit becomes too high for the 
czar. The gap between 0.7 deg/day and 1.2 deg/day leaves a 'donut hole' in the velocity 
space where pairwise tracklets will not be found. Recent improvements in our ability to 
screen false detections will allow this hole to be removed when Pan-STARRSl is observing 
away from the Galactic Plane. 

In the fast pairwise tracklet regime MOPS makes use of morphological information 
{e.g. moments) provided by the IPP as a proxy for trailing information. Quality cuts 
require similar morphology between the two detections and that the position angles (PA) 
calculated from the detection moments be aligned with the PA of the tracklet. 
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5.2. Detections 

The false detection rate delivered by the Pan-STARRSl IPP to MOPS is shown in 
Fig. [151 The rate is dominated by systematic image artifacts rather than statistical noise. 
At high galactic latitudes the transient detection rate is about the expected value from 
statistical fluctuations in the background but the rate increases dramatically as the field 
center approaches the galactic plane where the false detection rate is 10 — 50 x higher than 
expected. The increase in false detections with proximity to the galactic plane is due to 
image subtraction issues in regions with high stellar sky plane density and residual charge 
in the CCDs from bright stars. The transient rate in the wpi filter does not have data for 
galactic latitudes < 30 deg because this filter was solely used for solar system surveying 
and detection of moving objects. The solar system team does not survey into the galactic 
plane because we learned early-on that the false detection and tracklet rates are too high 
to efficiently detect solar system objects. 

Overall, the rate of transient detections at >5a is ~ 8200/deg^ — orders of magnitude 
higher than expected with Poisson statistics and with far more systematic structure, 
e.g. from internal reflections and other artifacts. This structure is especially problematic 
because, unlike Poisson-distributed false sources that are evenly distributed across the 
detector, systematic false sources tend to be clumped spatially and form many false tracklets 
that can clog and contaminate subsequent data processing. A large number of reported 
false detections for which the S/N is much lower than expected at a given magnitude could 
easily be removed with a fleld-dependent cut on S/N vs. magnitude (see Fig. [T6|l . 



Figure [T7] shows a representative sample of the types of false detections provided to 
MOPS by the IPP. Ac curate machine clas sification of these artifacts using methods such 



as those described by ( iDonalek et al. 



20081 ) can greatly reduce contamination and is under 



exploration by the PSISC. We use imaginative monikers like 'arrowheads', 'chocolate chip 
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cookies', 'feathers', 'frisbees', 'pianos' and 'UFOs' to classify them when iterating with the 
IPP on improving the detection stream. Many of the false detections are easily explained 
as internal reflections, ghosts, or other well-understood image artifacts but some are as 
yet unexplained. MOPS screens false detections using their morphological parameters but 
many survive the cuts into tracklet formation because they must be conservative so that 
diffuse comet detections are not removed. 

Finally, when the density of false detections spikes to rates that are difficult to manage, 
e.g. near a very bright star, we invoke a last-ditch spatial density filter to reduce the local 
sky-plane density to no more than 10,000 detections/deg^ in a circle of 0.01° radius. The 
filter eliminates false detections by discarding those with the highest PSF-weighted masked 
fraction until the detection density is manageable. 

Even real detections are often not well formed PSFs as illustrated in Fig. [181 Trailed 
detections often intersect with cell and chip gaps (a), PSF-like detections of slower moving 
objects have overlapping PSFs (2012 PF12), or sit on the boundary of a masked region 
{e.g. e). Whenever a detection is poorly formed it creates difficulties in astrometry and 
photometry and therefore in linking the detections into tracklets and tracks. 

Despite these problems Pan-STARRSl astrometry and photometry reported to the 
MPC for slow moving objects is currently excellent by contemporary st andards. The 
average RMS astrometric uncertainty is about 0.13" (iMilani et al.ll2012l ). The intrinsic 
photometric uncertainty in the calibratecj^ Pan-STARRSl dat a is < 10 mmag in t he (^pi. 



2OI2I). The 



rpi, and i-pi filters and ~ 10 mmag in the z-pi and y-pi filters (iSchlafiy et al. 
asteroid photometry is not yet as good as in the 5 primary Pan-STARRSl filters because 
1) much of the data is in the still to be fully calibrated w-pi filter and 2) the detections 



^Pan-STARRSl began reporting calibrated magnitudes in May 2012. 
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are often trailed enough to cause problems with the photometric fit. We expect to solve 
these problems soon and begin reporting wpi band photometry for moving objects with a 
calibration accuracy of a few tens of mmag. 



5.3. Tracklets 



Kubica et al.l ( 120071 ) provide detailed information on our algorithms and their 
performance for linking detections of moving objects on a single night into 'tracklets' or 
linking tracklets between nights into 'tracks'. Their techniques rely on the use of 'kd-tree' 
structures to provide fast and scalable performance. 

Their Table 1 shows that the algorithms yield nearly 100% tracklet creation efficiency 
for the S3M in the presence of the expected but random false detection rate of 250 deg~^. 
About 10-15% of the generated tracklets are MIXED or BAD but this tracklet accuracy is 
perfectly acceptable because false tracklets should not link together across nights to form 
a set of detections with a good orbit. The mixed tracklets cause no loss of objects because 
the tracklet creation is non-exclusive i.e. the same detection can appear in the correct 



Kubica et al 



( 120071 ) note that the technique 



CLEAN tracklet and in multiple other tracklets. 
still worked with high efficiency even without using the detections' orientations and lengths 
but at the expense of a much higher false tracklet rate. Furthermore, they showed that the 
system maintains its integrity even at > lOx higher random false detection rates. 

The algorithm for linking tracklets between nights uses a kd-tree to quickly identify 
candidate tracks. Each track is then fit to a quadratic in RA and Dec with respect to time 
to eliminate tracks based on their fit residuals. The number of candidate tracks increases 



exponentially with the number of detections as shown in Table 3 in iKubica et al. 



(l2007h 



while 98.4% of the correct tracks are properly identified only 0.3% of the tracks that pass 
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all the cuts are real. The false tracks are eliminated in the next step through initial orbit 
determination (lOD) described in the following sub-section. 

As discussed in the previous section, tracklet formation must take place in the presence 
of false detections due to statistical fluctuations on the sky background but also the 
usually much more numerous and problematic systematic false detections. We note that 



Kubica et al. 



(120071 ) performed their tests using 2-detection tracklets and with assumptions 
on the astrometric and photometric performance expected for the cosmetically clean images 
of an idealized Pan-STARRS4 survey. Their simulations are thus not directly applicable 
to the real Pan-STARRSl system that includes more noise and degraded astrometry. The 
Pan-STARRSl false detection rate is 10 — 50 x greater than expected at > 5(t and almost 
all of them are systematic. They pose problems for MOPS in that the detections "line up" 
well enough to pass the tracklet formation cuts and the tracklets represent objects with 
asteroid-like rates of motion. The high false tracklet rate significantly slows down derived 
object processing because of the resulting tremendous false track rate. 

The realized tracklet creation efficiency and accuracy was measured with synthetic 
detections injected into the real Pan-STARRSl data. In this data stream typically 4 
images are acquired at the same boresight with roughly a TTI separating each image in 
the sequence. We form tracklets of 2, 3, or 4 detections with some limits on the effective 
rate of motion depending on the number of detections in the tracklet. Table [8] shows that 
the tracklet creation efficiency is 99.98% — the algorithm correctly identifies almost every 
possible set of detections in a tracklet even in the presence of real detections and noise. 

The algorithm returns additional MIXED and BAD tracklets though at a small fraction 
of the CLEAN rate (see §3.11 and §3.111 for definitions). Tracklets are also identified for the 
non-synthetic detections i.e. detections of real objects or false detections. For the purpose 
of the efficiency and accuracy calculation both types of non-synthetic detections are 'noise' 
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as viewed from the perspective of the synthetic solar system objects. Thus, non-synthetic 
tracklets may be sets of detections of real solar system objects — but the numbers in Table [S] 
indicate that there is a problem. There are roughly 7x more non-synthetic detections than 
synthetic but, if our simulation and synthetic population are an accurate representation of 
reality, and all the non-synthetic detections are real, the two values should be roughly the 
same. 

Tracklets can only be formed from detections in images acquired within a reasonably 
short time frame with sequential detections of the same object separated by about a 
TTI. Linking detections into tracklets across longer time intervals creates a combinatoric 
explosion that can not be addressed by the current linking algorithms. On the other hand, 
it often happens that multiple tracklets for the same object are created on the same night 
separated by time intervals large compared to a TTI. The most common cause is overlap 
between adjacent fields acquired at times separated ^TTI. Even more problematic are 
'deep-drilling' situations where ^ 4 images are acquired at the same boresight in rapid 
succession. In this situation a real object moving between the images is not necessarily 
detected in every image due to e.g. chip and cell gaps, passing over a star, masking, 
photometric fluctuations, etc., and detections from the same object may appear in different 
tracklets. Furthermore, the detections within the tracklets may overlap in time due to 
astrometric fluctuations, {e.g. detections from images 1, 3 and 6 appear in tracklet A and 
detections from images 2, 5, and 8 appear in tracklet B with no detections in images 4 and 
7.) Thus, when there are > 4 exposures in a sequence at the same boresight we post-process 
the tracklet list with collapseTracklets. 

collapseTracklets merges co-linear tracklets within a single night or data set using a 
method similar to a Hough transform. The approximate sky plane location of each tracklet 
is determined at the mid-time of all tracklets in the set assuming that their motion is 
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linear in RA and Dec during the time frame (d and 5 respectively). Co-linear tracklets 
corresponding to the same object should have similar positions and motion vectors making 
them straightforward to identify as 'clumps' in {a, 6, a, 5)-space using a series of range 
searches implemented with a 4-dimensional kd-tree. Given candidate groupings of tracklets 
in the transform space, collapseTracklets attemps to merge them using straightforward 
RMS residual acceptance criteria, where there is a basic tradeoff between creating false 
merges vs. allowing too many duplicates. The cost of duplicate tracklets in the deep-drilling 
sequences is that if there are duplicate tracklets for the same object, a potential derived 
object will rightly use the distinct tracklets to create two distinct derived objects for the 
same object, and derived object discordance rejection (see §5.51) will disallow the derived 
objects on the grounds that multiple tracklets for an object cannot exist in the same set of 
fields. 

We adopted collapseTracklets essentially unmodified upon delivery from the 
development team at LSST, and beyond simple tests that verify basic capability, we have 
not attempted to maximize its efficiency with Pan-STARRSl deep-drilling cadences or data 
quality. Informal evaluations against the Pan-STARRSl 8-exposure Medium Deep (MD) 
sequences show a duplicate ratio of ~25% with essentially zero lost tracklets. 

In the time span from February 2011 through May 2012 MOPS created 1,513 x 10'^ 
tracklets of which 534 x 10^ were real (35%) and 345 x 10'^ were automatically attributed to 
numbered or multi-opposition objects (see §5.91) . The situation is much worse for tracklets 
containing just 2 detections: 11,691 x 10^ created with just 780 x 10^ attributed (6.7%). 
Assuming that the attributed:real ratio is the same as for the 3- and 4-detection tracklets it 
implies that only ~10% of 2 detection tracklets are real. 

In theory, the Pan-STARRSl-|-MOPS system is capable of detecting the highest proper 
motion stars within a single lunation. Consider the unusual case of Barnard's star with 
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,1 



an annual proper motioqll of ~ 10.3"/year or about 0.0282^yday . Pan-STARRSl's average 



2OI2I ) so that Barnard's star 



RMS astrometric uncertainty is about 0.13" flMilani et al. 
shows noticeable astrometric motion in about 5 days. Barnard's star would be provided 
as a transient detection to MOPS if the Pan- STARRSl IPP creat ed difference images 



20081). On each night MOPS 



using a static sky created from earlier images (iMagnier et al.l 
would create a 'stationary' transient and if the survey provided more detections 5 and 10 
days later MOPS would link the tracklets together. Thus, even in normal operations it 
is theoretically possible that MOPS could provide detections of high proper motion stars 
light years from the Sun. It takes little imagination to realize that MOPS could process 
'stationary' tracklets identified in three successive months to discover proper motion stars 
to even larger distances. The false detection rates for image differencing using a static sky 
should be no worse than pairwise image differences at the same S/N. 



5.4. Tracks 



Tracklets are linked into 'tracks' across previous nights over a time w indow of typic ally 



7-14 d ays using another kd-tree implementation called linkTracklets ( iKubica et al 



( 120071 )). The MOPS runtime configuration specifies how many days prior to the current 
night to search and how many tracklets are required to form a linkage. The tracklets search 
locates combinations of tracklets that collectively fit quadratic sky-plane motion to within 
some configured error. Of course actual asteroid motions are not truly quadratic, but within 
the time interval of track creation this approximation holds for most objects. 

Column 4 in table |9] shows that the tracklet linking efficiency is ~ 80% for most classes 
of solar system objects ranging from the inner solar system to beyond Neptune. The 



if 



http : //simbad . u-strasbg . f r/simbad/siiii-id?Ident=V[i'+V2500+0ph 
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linking efficiency is reduced from our grid simulations for several reasons: a) linkages are 
contaminated with false tracklets and fail orbit determination; b) Medium Deep sequences 
generate multiple tracklets for the same object, causing rejections in the discordance checks 
(see I5.5p : and c) we simply have not yet tuned linkTracklets' operational parameters to 
Pan-STARRSl's current survey mode and astrometric perfor mance. Early MO PS tests 



with synthetic and real data obtained from Spacewatch {e.g. iLarsen et al.ll200ll ) showed 



that the derived object efficiency can be increased to nearly 100% with a suitable survey 
strategy, false detection rate, and set of linkTracklets configuration parameters 

Table M does not show linkTracklets's accuracy but it is similar to the ^1% accuracy 



reported by 



Kubica et al. 



(120071 ) ■ i.e. >99% of all the tracks are not real. However, very 
few of the tracks survive the subsequent motion, acceleration, astrometric and photometric 
residual cuts, etc. Those that do pass all the cuts are subject to initial and differential orbit 
determination to eliminate essentially all the bad tracks as described in the following two 
subsections. 



5.5. Orbit Determination 

MOPS identifies tracks that represent a real or synthetic asteroid through 'orbit 
determination'. A six-parameter orbit is calculated for each track and only those that do 
not exceed a RMS residual requirement are accepted and passed on to the next processing 
phase. Orbit determination is a two-stage process within MOPS beginning with initial 
orbit determination (lOD) and followed by a least-squares differentially corrected orbit 
determination. Tracks that do not satisfy the residual requirement are discarded and their 
tracklets made available for use in other linkages. 

Occasionally there may be cases where a tracklet is included in multiple distinct 
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linkages whose fitted orbits produce RMS residuals that meet MOPS acceptance thresholds. 
When this occurs the tracklet cannot logically be included in both derived objects, so 
MOPS includes code 'discordance identification' software that identifies this situation and 
either 

• chooses one linkage as 'correct' if its RMS residual is significantly smaller than all 
others, or 

• rejects all linkages as bad since no correct linkage can be determined. 

Similarly, two distinct tracklets may share a common detection and be included in 
otherwise separate linkages with acceptable RMS residuals. In this case it is logically 
impossible for the detection to belong to two different objects, so the discordance 
identification routines select one linkage as described above or rejects them all if none 
stands out as correct. Tracklets in rejected linkages can be recovered later in processing 
through precovery (see §5.10p . 



5.6. Initial Orbit Determination (lOD) 

Orbit determination using a small number of observations over a short time interval 
has been studied for centuries, since Kepler first described the laws of planetary motion. 
Many techniques have been devised, each with its own advantages, and modern computers 
allow many different methods to be evaluated for a set of input detections. 

For MOPS's purposes, lOD is a module that produces six-parameter orbits given a set 
of detections. We irnpleme nted the 'OrbFit' package developed by the OrbFit Consortium 



( iMilani and Gronchi 



20101 ). Our selection criteria for an orbit determination package were 



1) efficiency at producing an orbit given a 'true' linkage; 2) orbit accuracy for observations 
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using Pan-STARRS cadences and astrometric uncertainties; 3) speed of computation; and 
4) availability and support. We tested the orbit determination software against millions of 
synthetic tracks to measure performance. As with other MOPS software, while we wanted 
high efficiency and accuracy in our orbit determination software, it was equally important 
to precisely measure the orbit determination efficiency. 

Table [9] shows that the lOD efficiency is essentially 100% for all classes of solar system 
objects, i.e. OrbFit successfully provides an orbit when provided a correctly linked set of 
tracklets (a track). 

5.7. Differential Correction 

MOPS attempts to improve the lOD by navigating the parameter space of orbit 
solutions to minimize the least-square RMS residual of the fitted detections. As with 
initial orbit determination, differential correction is a long-studied problem, spanning many 
decades of research, made routine by today's computer hardware. Modern differential 
correction techniques use precise models for the motion of the Earth and other large solar 
system bodies and are able to predict positions on the sky to within hundreths of an 
arcsecond many years into the future given enough input observations. 

Through a memorandum of understanding (MOU) with the Jet Propulsion Laboratory 
(JPL) Pan-STARRS obtained permission to use a subset of JPL's Solar System Dynamics 
(SSD) software to compute asteroid ephemerides and differentially corrected orbits. SSD is 
a workhorse of solar system analysis, used for guiding spacecraft missions to their asteroid 
and comet destinations and for assessing impact probabilities for high profile NEOs. 

Table [9] shows that the differential correction efficiency for objects that have passed 
lOD is essentially 100% for all classes of solar system objects. 
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Despite the overall effectiveness of the OrbFit lOD+JPL differential corrector 
combination, certain geometric configurations can lead to convergence to an incorrect orbit 
— the motion of objects in the small-solar-elongation sweetspots leads to dual solutions 
for orbits computed with short arcs. For short windows around the orbit computation, 
typically up to 30 days, the ephemeris uncertainty using the 'wrong' orbit is small enough 
that an attribution or precovery tracklet can still be found. With the additional arc length 
from the attribution or precovery the new orbit will collapse to a single solution. 



5.8. Derived Objects 

The tracklets and orbital parameters of tracks that survive the RMS residual cuts for 
orbit determination are stored in the MOPS database as a 'derived object'. A derived 
object represents a moving object 'discovered' by MOPS; in other words, the detections are 
believed to belong to the same body with an orbit that allows the body's motion to be 
predicted well enough to recover the body in an adjacent lunation to the discovery lunation. 

The derived object efficiency and accuracy is provided in table [91 The realized 
efficiency for the Pan-STARRSl survey is currently in the 70-90% range for most classes of 
solar system objects. Due to the high false detection rates the survey has con centrated on 



identi 



1991 



ying candida te NEO tracklets by their anomalous rates of motion {e.g. 



Rabinowitz 



Jedicke 



19961 ) and has temporarily abandoned the idea of creating derived objects. 
Despite a survey pattern, cadence and un-optimized MOPS configuration that are not 
particularly well-suited to creating derived objects, the system still achieves 70-90% 
efficiency for objects that appear in multiple tracklets within a lunation. 
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5.9. Attribution of known objects 

MOPS was designed to operate agnostically on all tracklets regardless of whether they 
correspond to known or unknown objects, i.e. no a priori information about known objects 
is used when creating MOPS tracklets, tracks or derived objects. The motivation was to 
create a final set of 'derived objects' from Pan-STARRSl-only data with good controls on 
observational selection effects. 

On the other hand, for the purpose of extracting more science from the data, and as 
an alternate means of measuring the system's detection efficiency and accur acy, we also 



attribute tracklets to known numbered and multi-opposition asteroids using iMilani et al. 



(l2012l )'s known_server module. They show that it has essentially 100% efficiency and 
100% accuracy using Pan-STARRSl data for those classes of objects and subsection §5.131 
discusses the realized Pan-STARRSl +MOPS detection efficiency characteristics using the 
known objects. 

As of October 2012 we have reported detections for approximately 240,000 numbered 
and 84,000 multi-opposition objects identified by known_server to the MPC These 
detections represent about 73% of the total 3.4 million detections, the remainder being 
mostly previously unknown asteroids. The fraction of known objects decreases with 
increasing y- magnitude such that 50% of the reported objects are unknown for V > 21.6 
(see Fig. [19]). 



5.10. Precovery & Attribution ('PANDA') of derived objects 



In keeping with the principle of MOPS agnosticism with respect to previously known 
objects we implemented the capability within MOPS to 'attribute' new tracklets each night 
to MOPS's derived objects. Recall that the time window for creating a MOPS derived 
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object is typically 7-14 days, but we want to associate individual tracklets from new data 
with existing derived objects if possible (attribution), and search the MOPS database for 
individual tracklets observed previously that did not form a derived object because not 
enough tracklets were observed at the time (precovery). Early Pan-STARRS work showed 
that with sufficiently high-quality astrometry, an orbit with arc length of typically 10-14 
days would have a prediction uncertainty small enough to locate the object in an adjacent 
lunation. Then, after a successful precovery or attribution extends the arc to beyond 30 
days, the orbit is secure. 



Milani et al. 



fl2012h 



The PANDA algorithm is a simplified version of that described by 
— in essence we integrate every derived object's motion to the time of observation of each 
image and compare the predicted location and velocity to all nearby tracklets to see if 
there is a match. If there is a match we attempt a differential orbit computation and if 
the resultant fit and residuals are within acceptable bounds we attribute the tracklet to 
the derived object. In practice, we pre-determine the locations and velocities of all derived 
objects at the begin ning, middle and end of the night and then use f ieldProximity 



(IKubica et al. 



20051 ) which interpolates their locations and velocities to each image time 
and uses a kd-tree implementation to make candidate associations between derived objects 
and tracklets. 

Table [10] shows that the attribution efficiency is ~93% but it is important to keep two 
points in mind: 1) this is the cumulative efficiency for all possible attributions on that night 
and 2) the statistics are dominated by the main belt asteroids. In regard to the first point, 
our efficiency determination software knows when a new synthetic tracklet is detected for 
an existing synthetic derived object but it does not account for the accuracy of the derived 
orbit's ephemerides at the time of observation e.g. the derived object's arc-length. It may 
be that the available attribution was for a derived object detected on each of the last 3 
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nights or it may be a derived object whose last observation was 3 years beforehand. Thus, 
we expect that attribution efficiency will decrease as a function of time. Even though 
the results in table [TD] are dominated by the main belt the attribution efficiency mostly 
increases with the semi-major axis of the object — it is relatively easy to attribute distant, 
slow moving objects because their sky plane density is low and their motion is mostly along 
a great circle, and it is difficult to attribute nearby objects for the opposite reasons. 

The mixed and bad attributions in table [10] are of particular concern (see §3.111 for 
definitions of the terms). In these types of attributions unassociated detections are added 
to the derived object, i.e. the derived object is being contaminated by noise or detections 
of other objects. Since the synthetic MOPS objects are intended to mimic the behavior 
of the real objects we assume that the real objects suffer the same 'contamination' levels 
(from real and synthetic detections and from false detections). Some level of contamination 
is always unavoidable but it will degrade the quality of the derived orbits so that future 
attributions are less likely to be real. 

There are several ways to mitigate the contamination including 1) tighter controls on 
the residuals of candidate detections added to derived objects or 2) post-facto 'scrubbing' 
of all derived orbits to identify outlying detections associated with derived objects. We 
never implemented these techniques and have essentially abandoned improvement of the 
attribution algorithm because it was designed for a survey pattern that would provide > 3 
tracklets for most observed asteroids in each lunation. 

The complement to attribution is 'precovery' in which historical unattributed tracklets 
are linked to new derived objects. Our precovery algorithm is essentially identical to 
attribution but applied backwards in time. The only enhancement is that when a successful 
precovery enhances the derived orbital element accuracy with a consequent reduction in 
the ephemeris errors we iterate on the precovery attempts until all possible precoverable 
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tracklets are associated with the new derived objectjfl We also allow for the precovery of 
tracklets intermediate in time between existing tracklets in the derived object. 

When the observed arc length of a derived object exceeds a configurable time span such 
that the prediction uncertainty over the entire survey is below some threshold, the object 
has reached a stage where the orbit solution is stable, all possible detections have already 
been associated with the object, and it no longer needs to be precovered. This optimization 
reduces the processing time required for precovery (see Fig. [20]) because over a long enough 
survey many objects will be 'retired' out of precovery when their orbits become sufficiently 
accurate. For synthetic objects, MOPS records all successful and failed precoveries so the 
pipeline can produce statistics on objects that have reached orbit stability but still have 
recoverable tracklets (see table ITT!). 

Fig. |2T] illustrates the sequence of events within MOPS for a hypothetical asteroid 
that is observed over three successive lunations. For synthetic objects, the 'paper trail' 
allows us to interrogate derived object performance and understand how and what kinds of 
objects get lost in the system. For each derived object 'event' (derivation, attirubtion, or 
precovery), the object's orbit is modified (hopefully improved) and recorded, allowing easy 
inspection of the orbit history for the derived object within MOPS. 

5.11. Orbit Identification 

When a new derived object {A') is created it is automatically checked for possible 
precoveries as described above. If the arc length is still short even after precovery it is 
possible that much earlier detections of the object A' appear in the database that can 

^The precovery algorithm is invoked any time any orbit changes, i. e. the attribution of a 
new tracklet to an existing derived objects also triggers the precovery algorithm. 
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not be precovered because the ephemeris uncertainty becomes too large when the current 
derived orbit is integrated backwards in time. On the other hand, those much earher 
detections may also have been incorporated into a derived object (A) that similarly can 
not be attributed to the current detections because integrating derived object A forward in 
time generates too large an ephemeris uncertainty. Thus, the same object may exist in two 
separate derived orbits that are unlinkable by our precovery and attribution techniques. 

Rigorous techniques exist for 'orbit identification' — associating short arcs of detections 
withi n an apparition to the same object observed wit h a short arc in another apparition 



.e.(/. 



Granvik and Muinonen 



2008 



Milani et al. 



20051 ) — but given the high astrometric 



quality we could assume for Pan-STARRS4, we adopted a simplified strategy that searches 
for neighboring orbits in orbit parameter space. 

If each of the two derived orbits are themselves reasonably accurate and precise then 
it is possible to identify the two derived objects as being identical simply by checking that 
the derived orbit elements are similar. Given that we expected to find > 10^ objects with 
Pan-STARRS and thousands of new derived objects each night we implemented the search 
for similar derived orbits in a kd-tree (orbitProximity). For any new derived object A' 
orbitProximity identifies a set of candidate derived orbits {X = {A, B, C, D}) for which all 
six orbit elements match A' to within the tolerances shown in table [T^l Then we attempt a 
differential correction ( §5.7p to an orbit including all detections from pairwise combinations 
of A' with each of the candidate orbits, e.g. {A, A'}, {B,A'}. If the differentially corrected 
orbit meets our RMS residual requirements then the two derived objects are merged into 
one derived object {e.g. A and A'). 

Our MOPS simulations for a Pan-STARRS4 system showed that the orbit identification 
efficiency and accuracy were close to 100% using the tolerances in table [T21 The realized 
orbit identification efficiency for Pan-STARRSl is ~26% as show in table [131 The 
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reduced efficiency is because the Pan-STARRSl cadence and astrometry yields orbits with 
uncertainties larger than our configured orbit identification tolerances. 

5.12. Data Rates & MOPS Timing 

Figure [20] shows execution times for various MOPS processing stages over one year 
of a 2-year MOPS simulation with the S3M. Processing times for per-night stages within 
the pipeline (generation of synthetics and tracklets) are generally constant, while stages 
that operate on the simulation's internal catalog of derived objects grows linearly or 
worse, depending on gross algorithmic considerations. Naively, the precovery stage exhibits 
quadratic growth, as the number of derived objects seen by MOPS and the number 
of precovery images to search both grow linearly as a function of simulation time, but 
optimizations in the precovery algorithm that limit the search windows reduces this growth 
below quadratic 0{n ■ log n)). 

The bimodality of timing results in Fig. [20] is due to the survey containing nights with 
either sweetspots (168 exposures) or opposition regions (660 exposures) or both. The slow 
growth of the synthetic generation stage is due to the increased integration time required 
to produce a position for an asteroid at the observation epoch from the survey start date. 
Recent MOPS versions eliminate this overhead by periodically propagating the epoch 
of S3M orbits to the current end-of-survey. Inter-night linking and orbit determination 
computation time grow slowly as the size of the MOPS database increases. 

5.13. Detection efficiency 

Figure [22] illustrates the realized tracklet detection efficiency of the combined 
Pan-STARRSl+MOPS systems as measured using asteroids attributed (or not) 
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by known_server (see §5.9p . The data in each filter were fit to the function e = 

-1 

where eo represents the maximum efficiency for bright but 
unsaturated detections, V and L are the apparent and hmiting V^-band magnitudes 
respectively, and w is the 'width' of the transition from maximum to zero efficiency. The 
limiting magnitude L is the magnitude at which the efficiency is 50% of the maximum 
value. The magnitudes in all bands were converted to V using a mean C+S class asteroid 
type as shown in table [H The function is a good representation of the efficiency near the 
limiting magnitudes and the scatter at bright values is simply due to low statistics. 

The maximum efficiency in all 4 filters in Fig. [22] is ~78% which is almost entirely 
driven by the effective Pan-STARRSl camera fill factor of ~80%. This implies that when a 
asteroid is imaged on a live and unmasked camera pixel it is detected with ~97% efficiency. 
The limiting magnitudes in ^fpi, rpi, and ipi are roughly equal at ~ 20.5 and any 
difference between the bands is due to unbalanced exposure times. The wide-band filter 
wpi goes ~ 1 \^-mag deeper through a combination of the 45 second exposure time and its 
~ 3x higher bandwidth. 

Figures [23] and [241 show that Pan-STARRSl+MOPS performance was roughly constant 
over the first year of operations. The dip in efficiency in the range 55700 <MJD<55800 is 
due to overly aggressive filtering of false detections during that time period. The fit to the 
limiting magnitude in each passband as a function of time is consistent with there being no 
change in any filter. 

In Fig. l25] we compare the predicted number of tracklets using the S3M to 
Pan-STARRSl tracklets reported to the MPC during the period 13 Aug 2011 through 11 
Oct 2011. For bright main-belt asteroids with wpi < 19, Pan-STARRSl reported about 
1/3 the predicted number. Fainter than wpi = 19 the reported rate drops further until the 
sensitivity limit is reached. Many factors account for the discrepancy between the predicted 



1 + exp ([V - L]/w 
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and realized rates: 

1. objects can be lost (primarily MBOs) in image differencing when observing away from 
opposition because of stationary points {i.e. the differencing subtracts all or part of 
the detection) 

2. tracklets submitted to the MPC might not yet be assigned a designation (a 'one night 
stand' tracklet) 

3. per-exposure live detector fraction may fall below the nominal 75% due to e.g. 
aggressive masking in regions with high stellar sky-plane density 

4. true night-to-night sensitivity in wpi is poorer than the constant model used for 
synthetics 

5. the synthetic S/N model overestimates faint-end performance 

6. the S3M may overstate the number of objects in some sub-populations. 

The situation is better for the NEOs that show better agreement between the model 
and Pan-STARRSl down to the sensitivity limit. The effect of (1) almost disappears 
because NEOs will rarely have stationary sky-plane motion, and (2) is reduced because of 
the NEO confirmation process for one-night tracklets coordinated by the MPC. 

6. Current Pan-STARRSl Performance 

As of late 2012, great strides have been made image processing and telescope scheduling 
to increase the yield of asteroids from the Pan-STARRSl survey. Better measurement and 
modeling of detector readout noise has resulted in greater detection sensitivity in image 
processing. Finer characterization of detection morphology by the IPP allows MOPS to 
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reject many classes of false detections upon ingest so that they do not contaminate the 
NEO tr acklet review pr ocess. The Pan-STARRSl Modified Design Reference Mission 



survey (jChambers 



2OI2I ) has increased the fraction of Sir survey time spent observing in 
quads. With the recent completion of the 37r static sky processing, there will be additional 
improvement in the system limiting magnitude of at least 0.4 mag as IPP moves from 
pairwise difference imaging to static sky differencing. 

Fig. [26] shows NEO discovery and submitted object totals for Pan-STARRSl and 
other NEO surveys. Although monthly totals can vary with weather losses, since late 2010 
the Pan-STARRSl NEO discovery rate has been increasing. Using IPP's morphological 
characterization and a MOPS comet candidate review procedure similar to that for NEOs, 
Pan-STARRSl has become a capable comet finder, discovering 30 to date and 8 in October 
2012 alone. 



7. Availability & Ongoing Development 

The MOPS software is available under the GNU General Public License Ver- 
sion 2, and can be retrieved by sending a query to the PSl Science Consortium 
(http://www.pslsc.org). MOPS has become a large and somewhat unwieldly package. 



employing (too many) different programming languages preferred by third-party software 
(not all available directly from the PSISC), and containing many installation dependencies 
for Perl and Python modules. Fortunately, documentation and development notes are 
available via the PSISC web site. 

Not all MOPS subcomponents are freely available and some are not integrated into 
the MOPS software distribution. In particular, the JPL Solar System Dynamics package 
that performs MOPS differential correction must be obtained directly from JPL. Other 
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packages such as OrbFit and the SLALIB positional astronomy hbrary must be downloaded 
separately. 

Approximately 15-20 full-time equivalent (FTE) years have gone into the development 



of MOPS and the S3M, and development continues in supp ort of Pan-S 



NEO survey and upcoming NEO surveys such as ATLAS (ITonry et al. 



ARR Sl's ongoing 



2OI2I). Prior 



to h andling 



e.g. 



Gehrels 



ive Pan-STARRSl data, MOPS processed testing data from Spacewatch 



1991 



Gehrels and Jedicke 



Asteroid Light Curve Survey (TALCS; 



19961) and raw dat a in support of the Thousand 



Masiero et al. 



20091). MOPS has additionally been 



instrumental in providing simulation results for Pan-STARRS, LSST and ATLAS, and as a 
research tool for many graduate students and postdocs. 



8. Future improvements 

The future of MOPS likely falls in several areas: 

• Optimizing execution speed for next-generation Pan-STARRS4-like data volumes 
through reduction of database I/O and greater in- memory processing. 

• Improved MOPS installation and configuration for non-Pan-STARRS surveys. 

• Modeling of photometric artifacts, e.g. partially masked PSFs and trails. 

• The linear processing model employed by MOPS is appropriate for a mature survey 
but can be inadequate when evaluating production parameters. It is difficult to merge 
separate MOPS databases created under different runtime configurations into a single 
coherent database. Future versions of MOPS should be more agile in its ability to 
perform non-linear processing and merge datasets. 
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• Increased fidelity of syntlietics. For example, although our simulations coarsely avoid 
the moon they do not incorporate the effects of varying sky sensitivity or effects due 
to the presence of bright objects in the field of view. In one sense, however, these 
effects are really an aspect of detector sensitivity and not MOPS itself, and they 
can be integrated via per-field detection detection efficiency parameters if they are 
measured by IPP. 

• Fill-factor modeling. There is experimental code to use the mask of Fig. E] as a crude 
live-pixel server (LPS, §3.121) to simulate detections lost on dead areas of the detector. 

• Simplified processing for large-scale simulations, e.g. two-body ephemerides instead of 
a perturbed dynamical model. 

9. Summary 

MOPS has proven to be a capable tool for detecting moving objects in the 
Pan-STARRSl transient detection stream. Despite the differences between the targeted 
Pan-STARRS4 capability and current Pan-STARRSl performance, we have deployed 
MOPS effectively in our search for NEOs and comets and in characterizing the main belt. 
Importantly, we have estimates of our efficiency of object detection so that we can provide 
a foundation for large-scale population studies. 

Additionally, we have made progress toward our design goal of creating a system 
that can detect objects and compute their orbits with > 99% efficiency for most solar 
system populations when provided next-generation survey astrometry and data quality. 
Our experience with Pan-STARRSl shows that the desired performance is challenging to 
achieve with real telescope data. 

With nearly one year remaining in the Pan-STARRSl survey, we expect further 
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broad improvements in astrometry, photometry, difference imagii ig and sensitiv ity in the 



2OI2I) is under 



final Pan-STARRSl data release. The Pan-STARRS2 telescope (iBurgettI 
construction adjacent to the Pan-STARRSl facility and will have a raw detector sensitivity 
and optical performance exceeding that of Pan-STARRSl. Operating modes currently 
under consideration for the two telescopes have them working in tandem, increasing 
performance even further. Combined with optimization and tuning of the Pan-STARRSl 
MOPS pipeline, we look forward to applying the full capability of MOPS to this data and 
to the subsequent riches that will lie in this catalog of moving objects. 
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Table 1. MOPS Pan-STARRSl Filter Transformations* 



Transform 


Solar 


Mean S+C 


texp (s) ^ 


V-gpi 


-0.217 


-0.28 


43 


V-rpi 


0.183 


0.23 


40 


V-ipi 


0.292 


0.39 


45 


V-zpi 


0.311 


0.37 


30 


V-ypi 


0.311 


0.36 


30 


V-wpi 


0.114 


0.16 


45 



Note. — Transformations from Johnson V to the six Pan-STARRS filters are provided for both 
a standard solar spectrum and asteroids with a mean S-l-C spectral type. These transformations 
have been implemented by the Minor Planet Center and AstDyS. 



Tonrv et al 



(|2012l ) 



^Exposure times as of October 2012. 
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Table 2. Pan-STARRSl MOPS Hardware Components (see Fig. [2]) 



Item 


Purpose 


Number 


Total Capacity 


Disk 


Database Storage 


2 X 10 TB 


20 TB^ 


Disk 


Administrative 


2 X 2.7 TB 


5.4 TB 


CPU 


Cluster processing 


8x4 cores 


32 cores (~100 GFLOPS) 


CPU 


Administrative'^ 


2x4 cores 


8 cores (~25 GFLOPS) 


Network Switch 


Network 


32 ports 


1 Gbps/port 



All database storage employs RAID6 for data integrity. 

^Administrative functions include a user console, pipeline and Condor workflow management, 
and the web interface. 
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Table 3. MOPS Database Storage Requirements 



Data Component PS4 Estimated (GB) PSl Estimated (GB) PSl Actual^ (GB) 



Fields 


0.2 


0.06 


0.02 


3(7 Detections 


500,000 


150,000 


N/A2 


5cr Detections 


500 


150 


122 


Derived Object Parameters 


1,000 


300 


0.1523 


Synthetic Object Parameters 


1,100 


330 


2.43^ 


Tracklets 


1,200 


360 


2.19^ 


Image Postage Stamps 


N/A 


N/A 


1,800^ 



^As of October 2012, 2.5 years into the 3.5-year Pan-STARRSl survey. 

^Processing of transients below 3a confidence is currently untenable for Pan-STARRSl MOPS 
due to the systematic false tracklet rate. 

^MOPS produces derived object parameters for only a small subset of Pan-STARRSl data until 
the pipeline can be tuned for Pan-STARRSl performance. 

^Pan-STARRSl operations uses a 1/10 sampled synthetic solar system model (S3M). A final 
Pan-STARRSl processing of its transient catalog will include a full S3M. 

^Pan-STARRSl single-exposure sensitivity reduces the actual tracklet count on top of the reduc- 
tion from a shorter survey than Pan-STARRS4. 

^Image postage stamps (200 x 200 pixels) are stored on a network file system, not in the MOPS 
database. They are locatable using database records. 



- 61 - 



Table 4. S3M Grid Model Semi-Major Axis Number Distribution (see Fig. [7]). 



Semi-major axis range (AU) 


Number of objects 


0.75-1.5 


50,000 


1.5-6.0 


200,000 


6.0-32 


50,000 


32-50 


200,000 


50-500 


50,000 


500-5000 


25,000 



Note. — The number of objects in each pseudo-logarithmic range and their absolute magnitudes 
were selected so that the sky-plane density of detected objects in the simulations will be similar. 
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Table 5. Tracklet (intra-night) Efficiency* & Accuracy* in the 2-year Pan-STARRS4 

MOPS full-density simulation. 



Avail. 


Clean % Unfound 


% 


Mixed 


% Bad 


% 


Non-syn. 


24056199 


24056191 100.0 8 


0.0 


417774 


1.7 553475 


2.3 


1854648 



Note. — This simulation used a realistic solar system model (S3M) and "next-generation" 
(0.01") astrometric uncertainty (about 10 x better than delivered by the best contemporary 
surveys like Pan-STARRSl). The tracklet terminology in the headings is described in detail 

in gam 

*See §3.11 and §3.111 for definitions of efficiency and accuracy. 
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Table 6. Derived Object Efficiency for eight^ different classes of synthetic solar system 
objects in the 2-year Pan-STARRS4 MOPS full-density simulation. 



Object Avail.* Clean* % Pass % Pass % 
Class^ (Linked) lOD^ Diff.^ 



NEO 


5203 


4994 


96.0 


4924 


94.6 


4924 


94.6 


MBO 


2043584 


2032676 


99.5 


2029618 


99.3 


2029618 


99.3 


TRO 


56214 


56061 


99.7 


55923 


99.5 


55923 


99.5 


CEN 


557 


556 


99.8 


556 


99.8 


556 


99.8 


JFC 


551 


546 


99.1 


541 


98.2 


541 


98.2 


LPC 


1714 


1713 


99.9 


1584 


92.4 


1584 


92.4 


SDO 


2289 


2286 


99.9 


2281 


99.7 


2281 


99.7 


TNO 


12719 


12710 


99.9 


12686 


99.7 


12686 


99.7 



Note. — This simulation used a realistic solar system model (S3M) and "next-generation" (0.01") 
astrometric uncertainty (about 10 x better than delivered by the best contemporary surveys like 
Pan-STARRSl). The tracklet terminology in the headings is described in detail in ^3. Ill 

*See ^3.11 and ^3.111 for definitions of efficiency and accuracy. 

^The impactor and hyperbolic models are omitted because they did not exist in the S3M at the 
time of this simulation. 

^NEO - Near Earth Objects; MBO - Main Beh Objects; TRO - Trojans; JFC - Jupiter Family 
Comets; LPC - Long Period Comets; CEN - Centaurs; SDO - Scattered Disk Objects; TNO - 
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Trans-Neptunian Objects 

^the number (and percentage) for which initial orbit determination was successful. 

^the number (and percentage) for which a differentially corrected orbit was successfully computed. 
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Table 7. Derived Object Efficiency for the 4-year Pan-STARRSl MOPS NEO-only 

simulation. 



Object 


Avail.* 


Clean* 


% 


Pass 


% Pass % 


Classi 




(Linked) 




lODi 


Diff.2 


NEO 


11515 


11025 


95.7 


9404 


81.7 9404 81.7 



Note. — This simulation used a realistic population of ~ 250, 000 NEOs with Pan-STARRSl-like 
astrometric uncertainty (0.1") simulation and a 100% fill-factor detector. This simulation was run 
during commissioning of the Pan-STARRSl telescope before the actual Pan-STARRSl survey was 
implemented. 

*See ^3.11 and ^3.111 for definitions of efficiency and accuracy. 

^the number (and percentage) for which initial orbit determination was successful. 

■^the number (and percentage) for which a differentially corrected orbit was successfully computed. 
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Table 8. Tracklet Efficiency* & Accuracy* in Pan-STARRSl Observing Cycle 143 

(2011-08-13 through 2011-09-10)"^ . 



Avail. 


Clean 


% 


Unfound 


% 


Mixed 


% 


Bad % 


Non-syn. 


246596 


246558 


100.0 


38 


0.0 


1969 


0.8 


8884 3.6 


1429073 



Note. — Measured performance for a single observing cycle of real Pan-STARRSl tele- 
scope pointings using realistic astrometric uncertainty (0.1") and a simplistic 75% fill-factor 
model. 

*See §3. II and §3.111 for definitions of efficiency and accuracy. 
^Excluding fields within 15° of the Galactic Plane. 
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Table 9. Cumulative derived object efficiency for 10 different classes of synthetic solar 
system objects in Pan-STARRSl Observing Cycle 143 (2011-08-13 through 2011-09-10). ^ 



Object Avail* Clean* % Pass % Pass % Dup> % 
Class^ (Linked) lOD^ Diff.^ 



IMP 


88 


68 


77.3 


67 


76.1 


67 


76.1 





0.0 


NEO 


24 


18 


75.0 


18 


75.0 


18 


75.0 





0.0 


MBO 


18301 


14608 


79.8 


14551 


79.5 


14461 


79.0 


6 


0.0 


TRO 


662 


538 


81.3 


538 


81.3 


522 


78.9 





0.0 


CEN^ 








0.0 





0.0 





0.0 





0.0 


JFC 


22 


18 


81.8 


18 


81.8 


18 


81.8 





0.0 


LPC 


50 


39 


78.0 


36 


72.0 


36 


72.0 





0.0 


SDO 


7 


1 


14.3 


1 


14.3 


1 


14.3 





0.0 


TNO 


17 


15 


88.2 


15 


88.2 


15 


88.2 





0.0 


HYP 


32 


27 


84.4 


19 


59.4 


19 


59.4 





0.0 



Note. — Measured performance for a single observing cycle of real Pan-STARRSl telescope 
pointings using realistic astrometric uncertainty (0.1") and a simplistic 75% fill-factor model. 

*See ^3.11 and ^3.111 for definitions of efficiency and accuracy. 

^Excluding fields within 15° of the galactic equator. 

^IMP - Earth Impactor; NEO - Near Earth Objects; MBO - Main Beh Objects; TRO - Trojans; 
JFC - Jupiter Family Comets; LPC - Long Period Comets; CEN - Centaurs; SDO - Scattered Disk 



Objects; TNO - Trans-Neptunian Objects; HYP - Hyperbolic (Interstellar) Objects 

^the number (and percentage) for which initial orbit determination was successful. 

^the number (and percentage) for which a differentially corrected orbit was successfully computed. 

^the number (and percentage) of duplicate derivations. The duplicates contain non-identical but 
intersecting sets of detections of the same object. 

^Centaurs were inadvertently omitted from the S3M for this study. 
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Table 10. Attribution Efficiency* & Accuracy* in Pan-STARRSl Observing Cycle 143 

(2011-08-13 through 2011-09-10). + 



Avail* 


Clean* % 


Mixed* % 


Bad* % 


Non-syn.* 


9684 


8971 92.6 


5 0.1 


26 0.3 


3625 



Note. — Measured performance for a single observing cycle of real Pan-STARRSl telescope 
pointings using realistic astrometric uncertainty (O.l") and a simplistic 75% fill-factor model. 

*See ^3.11 and ^3.111 for definitions of efficiency and accuracy. 

^Excluding fields within 15° of the Galactic Plane. 
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Table 11. MOPS Precovery Efficiency* & Accuracy* in Pan-STARRSl Observing Cycle 

143 (2011-08-13 through 2011-09-10). ^ 



Avail* 


Clean* % 


Mixed* % 


Bad* % 


Non-syn* 


8557 


6635 77.5 


14 0.2 


23 0.3 


1966 



Note. — Measured performance for a single observing cycle of real Pan-STARRSl telescope 
pointings using realistic astrometric uncertainty (O.l") and a simplistic 75% fill-factor model. 

*See ^3.11 and ^3.111 for definitions of efficiency and accuracy. 

^Excluding fields within 15° of the galactic equator. 
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Table 12. Orbit Identification Matching Tolerances 



Orbit 


Tolerance 


Element^ 




perihelion 


0.1 AU 


eccentricity 


0.05 


inclination 


0.1° 


ascending node 


1° 


arg. of perihelion 


1° 


time of perihelion^ 


0.1 



Note. — The maximum difference between two derived object's orbits before they are tested for 

consistency with being the same object. 

^Cometary form. 

^In terms of fractional orbit period. 
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Table 13. MOPS Orbit Identification Efficiency* & Accuracy* in Pan-STARRSl 
Observing Cycles 143 & 144 (2011-08-13 through 201 1-10-10). + 



Avail* 


Found* 


% Mixed* 


% 


Bad* % Non-syn. 


40 


7 


17.5 1 


2.5 


2 5.0 22 



Note. — Measured performance for a single observing cycle of real Pan-STARRSl telescope 
pointings using realistic astrometric uncertainty (O.l") and a simplistic 75% fill-factor model. 

*See ^3.11 and ^3.111 for definitions of efficiency and accuracy. 

^Excluding fields within 15° of the galactic equator. 
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Fig. 1. — High level flowchart for the Pan-STARRSl Moving Object Processing System 
(MOPS). The data processing proceeds from left to right in the flgure, then repeats for 
every additional night of data ingested by MOPS. 
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MOPS PS1 Cluster Organization 
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Fig. 2.— MOPS Pan-STARRSl cluster configuration. 
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MOPS Database Design 
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Fig. 3. — MOPS database schema. S3M columns indicate synthetic data for efficiency 
assessment. 
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Fig. 4. — Pan-STARRSl MOPS web interface showing the sky map for the night of MJD 
56210. On the left are hsts of observing blocks executed by the telescope for the night. 
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Synthetic PS1 Detector Mask 
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Fig. 6. — Synthetic 7.0 deg^ field-of-view static focal plane mask for Pan-STARRSl MOPS 
simulations that simulates losses due to chip gaps and other masked area. The horizontal 
and vertical bands simulate the chip gaps while the remaining squares represent losses due 
to bad cells or those used for guide stars. The squares do not correspond to actual masked 
area on the real Pan-STARRSl detector; they exist simply to simulate additional area loss. 
The actual distribution of masked area on the real Pan-STARRSl detector is very different 
— there is more power on small scales than shown here. The overall simulated camera fill 
factor is fixed at 75% in agreement with the measured values shown in Fig. |23l 
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MOPS Grid Tracklet Efficiency vs. Log(Semi-Major Axis a) 
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Fig. 7. — Tracklet creation efficiency as a function of semi-major axis for a one lunation 
MOPS simulation using S3M grid objects. Note the truncated Y-axis range of 0.95 to 1.0. 
Fifteen fast-moving objects with sky-plane rates that exceed the configured limits for the 
production Pan-STARRSl MOPS account for the reduced efficiency at small semi-major 
axis. These hmits will be tuned upon reprocessing the full Pan-STARRSl survey. 
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MOPS Grid Tracklet Efficiency vs. Eccentricity e 
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Fig. 8. — Tracklet creation efficiency as a function of eccentricity for a one lunation MOPS 
simulation using S3M grid objects. Note the truncated Y-axis range of 0.95 to 1.0. 
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Fig. 9. — Tracklet creation efficiency as a function of inclination for a one lunation MOPS 
simulation using S3M grid objects. Note the truncated Y-axis range of 0.95 to 1.0. 
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MOPS Grid Derived Object Efficiency vs. Log(Semi-IUIajor Axis a) 
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Fig. 10. — Derived object creation efficiency as a function of semi-major axis for a one 
lunation MOPS simulation using S3M grid objects. The total efficiency is 97.5%. Losses at 
small semi-major axis are due to conservative acceptance limits in the current Pan-STARRSl 
production MOPS that reject objects with large sky-plane accelerations and/or large RMS 
residuals. These limits will be tuned upon reprocessing the full Pan-STARRSl survey. 
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MOPS Grid Derived Object Efficiency vs. Eccentricity e 
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Fig. 11. — Derived object creation efficiency as a function of eccentricity for a one lunation 
MOPS simulation using S3M grid objects. The failures at high eccentricity (e > 0.99) occur 
in initial orbit determination (lOD) for grid objects with semi- major axis a < 5. These 
failures can be corrected by modifying lOD acceptance parameters. 
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Fig. 12. — Derived object creation efficiency as a function of inclination for a one lunation 
MOPS simulation using only grid objects. 
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MOPS Grid Tracklet Efficiency vs. Declination 
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Fig. 13. — Tracklet creation efficiency as a function of declination for a one lunation MOPS 
simulation using PSl telescope pointings and S3M grid objects with Pan-STARRSl realistic 
(0.1") astrometric uncertainty. The uneven declination coverage is due to uneven coverage of 
the Pan-STARRSl survey during the lunation (14 Aug 2011 through 12 Sep 2011). Several 
fast- moving objects in the declination=89° bin were lost due to overly conservative tracklet 
acceptance parameters. 
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MOPS Grid Derived Object Efficiency vs. Declination 



-20 



20 



40 



60 



n = 8983 

— I 

80 



n Unfound 
■ Found 




20 40 

Declination, deg 



Fig. 14. — Derived object creation efficiency as a function of decfination for a one lunation 
MOPS simulation using S3M grid objects with Pan-STARRSl realistic (0.1") astrometric 
uncertainty. The uneven declination coverage is due to uneven multi-night coverage of the 
Pan-STARRSl survey during this particular lunation (14 Aug 2011 through 12 Sep 2011). 
While the total grid derived object efficiency is 97.5%, we beheve this number can be im- 
proved further by tuning of MOPS linking parameters. 
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Fig. 15. — Number of transient detections reported to the MOPS by the IPP as a function 
of galactic latitude in all six Pan-STARRSl filters. Almost all the transient detections are 
false. 
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Pan-STARR51 iv-band Transient Detections, MJDSHiCl^ 
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Fig. 16. — Pan-STARRSl transient detections reported to MOPS by the IFF as a function 
of S/N and wpi magnitude on MJD 56028. Grey dots are detections in tracklets that are 
real asteroids. Black dots are all other detections. The 'fork' in the distribution at wpi ~ 18 
is due to a change in observing conditions that caused distortion of FSFs for bright objects. 
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Nn-STARRSI Systematic False Detection Gallery 
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Fig. 17. — A sample of false source detections delivered to MOPS. The cause of the detection 
is provided under the image when known. The source detections are in the positive image 
(white pixels). Each image is a 200 x 200 pixel difference image with the original source 
detection in the center, i.e. the difference between two successive images acquired at the 
same boresight. Solid dark grey regions represent gaps between CCDs or cells or masked 



regions. 
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Fig. 18. — A sample of real object detections delivered to MOPS via difference imaging. 
Each image is a 200 x 200 pixel difference image with the positive source detection in the 
center. Many detections lie near artifacts or in heavily masked regions or are difficult to 
distinguish from systematic false detections. Objects a through / are unknown asteroids. 
Solid dark grey regions represent gaps between CCDs or cells or masked regions. 
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Fig. 19. — Overlapping detected and known asteroid distributions for a single night of 
Pan-STARRSl solar system observing. 
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Fig. 20. — Processing times for different MOPS pipeline stages over one simulated year of a 
two-year simulation using Pan-STARRS4 data volumes. Some data prior to MJD 54850 was 
discarded because MOPS pipeline software was not separating processing times for synthetic 
and tracklet stages. 
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MOPS Derived Object Processing 
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Fig. 21. — MOPS derived object processing overview for a single hypothetical asteroid, 
showing observed tracklets and the created derived object. In lunation 132, no derivation is 
available because there is only a single tracklet observed. In lunation 133, three tracklets are 
observed in a 10 day window, so a derived object is created and the precovery in lunation 132 
becomes available. In lunation 134, two attributions are available using predictions from the 
object's derived orbit. Numbers in circles show the sequence of 'paper trail' records inserted 
into the MOPS database. 
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Pan-STARRS1 Estimated Moving Object Detection Sensitivity, Night 55803 
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Fig. 22. — Pan-STARRSl moving object detection efficiency on a single night (MJD 55803) 
for detections of known numbered and multi-opposition asteroids in each filter as a function of 
V magnitude. The data in each filter were fit to the function e = 1 + exp ^[V^ — L]/w^ 
as described in §5.131 The vertical dashed line is at the V magnitude where the efficiency 
drops to 50% of the maximum (L). The vertical dotted lines provide the range [L — w,L + w] 
over which the efficiency drops quickly. 
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Pan-STARRSI Estimated Moving Object Detection Efficiency, iUJDs 55649-56054 
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Fig. 23. — Pan-STARRSI moving object detection efficiency for brigfit non-saturated de- 
tections as a function of MJD corresponding to tlie time period from approximately Feb 
2011 through Jun 2012. Around MJD 55820, MOPS began employing less-aggressive false 
detection filtering of IPP transient detections, boosting per-exposure detection efficiency to 
75% consistently. The dashed line is a spline-fit to the data. 
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Pan-STARRS1 Estimated Moving Object Detection Sensitivity, MJDs 55649-56054 
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Fig. 24. — Pan-STARRSl nightly limiting V magnitude in each of the four main filters 

used to detect moving objects as a function of MJD corresponding to the time period from 
approximately Feb 2011 through Jun 2012. The dashed lines represent fits to the passband 
data as a function of time. 
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PS1 i/ir-band Simulated vs. Measured Performance 
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Fig. 25. — Comparison of MOPS synthetic tracklet counts to Pan-STARRSl submitted 
tracklets during the period 13 Aug 2011 through 11 Oct 2011 for (left) main belt objects 
and (right) near-Earth objects. 
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Pan-STARRS1 NEO Discoveries by lUlonth 
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Fig. 26.— Pan-STARRSl discovery statistics as of November 2012 (from the Minor Planet 
Center) . 
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